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Abstract. 

Knowledge of the binary population in stellar gr oupings provides important information about the outcome of the star forming 
process in different environments (e.g., Bla auw|[l99 ll and references therein). Binarity is also a key ingredient in stellar popu- 
lation studies, and is a prerequisite to calibrate the binary evolution channels. In this paper we present an overview of several 
commonly used methods to pair individual stars into binary systems, which we refer to as pairing functions. These pairing func- 
tions are frequently used by observers and computational astronomers, either for their mathematical convenience, or because 
they roughly describe the expected outcome of the star forming process. We discuss the consequences of each pairing function 
for the interpretation of observations and numerical simulations. The binary fraction and mass ratio distribution generally de- 
pend strongly on the selection of the range in primary spectral type in a sample. The mass ratio distribution and binary fraction 
derived from a binarity survey among a mass-limited sample of targets is thus not representative for the population as a whole. 
Neither theory nor observations indicate that random pairing of binary components from the mass distribution, the simplest 
pairing function, is realistic. It is more likely that companion stars are formed in a disk around a star, or that a pre-binary core 
fragments into two binary components. The results of our analysis are important for (i) the interpretation of the observed mass 
ratio distribution and binary fraction for a sample of stars, (ii) a range of possible initial condition algorithms for star cluster 
simulations, and (iii) how to discriminate between the different star formation scenarios. 
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1. Introduction 

Observations and si mulations suggest that most stars form in 
binary systems (e.g. Duquennov & Mayor 199 It iMason et alj 
199 81 IGoodwin & Kroupal 120051; IKobulnickv & Frverl [2007t 
Kouwenhoven et al.1 120051120071)1 IGoodwin et alj 120071) . and 
that a subst antial fraction are part of a triple or higher-order 
system (e.g. Tokovinin & Smekhov 2002; Correia et al. 2006; 
Tokovinin et all 120061: iHu et al]|2008l) . Multiplicity is thus a 
fundamental property of the star forming process. Detailed 
knowledge of a young binary population can be employed to 
study the outcome of star formation, and consequently the star 
formation process itself. 

Surveys for binarity have indicated that the proper- 
ties of the binary population are a function of the spec- 
tra l type of the p rimary star. Practic ally all O-type stars 



(Mas on et al.l ll998) and B/A-type stars dShatskv & Tokov inin 
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20021; IKobulnickv & Frvei] 120071; iKouwenhoven et ai]|2007b|) 
are found in binary or multiple systems. lAbt & Levvl ( 1 19761) 
report a multiplicity fraction of 55% among F3-G2 stars, 
and in their CORAV E L spe ctroscopic study of F7-G9 stars 
Duquennoy & Mayon (1199 11) find a binary fraction of ~ 
60%. The binary fraction among M-type stars is 30 



40% (Fischer & M arcvlll992l;lLeinert et alJl997tlReid & Gizis 



tion decreases to 10 - 20% (e.s.. Gizis et a 


l.ll2003tlClose et al. 


2003 


; Bouv et al. 


2003 


; Bureasser et al. 


2003; Siegleretal. 


2005 


; Ahmic et al. 


2007 


; MaxtedetalJ|2008). 



In this paper we discuss in detail several methods of pair- 
ing individual stars into binary stars. We refer to the latter algo- 
rithms as "pairing functions". Several of these have a physical 
motivation, others are discussed because of their mathemati- 
cal simplicity. All these pairing functions have in common that 
they are frequently used in literature. The main goal of this 
paper is to explain the consequences of adopting a particular 
pairing function when doing a numerical simulation, or when 
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interpreting observations. A good understanding of the conse- 
quences of each pairing algorithm for the binary population is 
important for 

- The interpretation of observations. Are the measured prop- 
erties (e.g. mass ratio distribution and binary fraction) for a 
certain sample representative for the population as a whole? 
What is the role of selection effects, and how can the differ- 
ent pairing functions be distinguished? 

- Initial conditions for sim ulations. In N-body s imulations, 



such as STAR LAB (e.g.. IPorteeies Zwart et alJl200lh and 
NBODY6 (e.g. jAarsethl 19991) . a mass ratio distribution in- 



dependent of primary mass is often adopted. What are the 
consequences of this approach? 
- Star formation. Which pairing functions are expected from 
the different star forming scenarios? Wh at numerical sim 



ulations of clustered star formation (e.g. jBate e t al. 2003 



Bate|[2008) predict? Do we expect random pairing from the 



initial mass function? Which observations are necessary to 
be able to distinguish between these scenarios? 

This paper is organised as follows. In §|2]we introduce the ter- 
minology used and discuss our assumptions. In §|3]we discuss 
the mass distribution, and in §|4]the mass ratio distribution. We 
briefly describe the origin of the pairing between binary stars 
in §|5] The major part of this paper is §|6] where we discuss 
the different binary populations resulting from the choice of 
the mass distribution and pairing function. The specific differ- 
ences are discussed in detail in §|7] and the dependence on the 
generating properties in §|8] Observational complications and 
a strategy to recover the pairing function are described in §|9] 
Finally, we summarise our results in § [T0l 

2. Method and terminology 

We study the differences of the various methods of pairing 
individual stars into binary systems by analysing numerically 
simulated binary populations. For each binary system, we re- 
fer to the most massive component M\ as the primary star, 
and the least massive component Mi as the companion star. 
Our adopted definition is purely based on the current mass of 
the components, i.e., irrespective of their relative luminosity or 
initial mass. We define the mass ratio as q = MijM\, so that 
< q < 1 . The total mass is denoted M T - M\ + M 2 for a 
binary system, and Mj = Ms for a single star with mass M$ . 
In several cases, we construct binary systems from star form- 
ing cores of mass Mq with star forming efficiency e, so that the 
total mass of the resulting objects is Mj = eMc- 

In our simulations each star is given a mass. A subset of 
the stars is assigned a companion, the other stars remain single 
stars. The companion is given a mass according to a pairing 
algorithm. We refer to the algorithm that is used to combine 
individual stars into binary systems as the pairing function. 
The pairing function of a binary population may for example 
be random pairing of both companions from the mass distri- 
bution. In the particular case of random pairing, primary and 
companion are swapped, if necessary, so that the primary is the 
most massive star. 



Depending on which pairing function is used, the mass of 
the primary star, and in several cases the companion star, is 
drawn from a mass distribution /m(M). We refer to fiwiM) 
as the generating mass distribution. We denote the resulting 
mass distributions for primary stars, companion stars, sys- 
tems, and single stars with fM,(M), fM 2 (M), fM T (M), and 
/m s (M), respectively. Note that the distributions over primary 
and companion mass in a stellar grouping are never inde- 
pendent, /m 1 ,m 2 (M 1 ,M 2 ) + /m 1 (Mi)/m 2 (M 2 ), as by definition 
M\ > M->_. The resulting mass distribution for all single and 
primary stars is denoted with fu^ s (M), and the mass distribu- 
tion for all individual objects (i.e., singles, primaries and com- 
panions), is denoted with / a ii(M). The distribution f d \\{M) that 
is present immediately after star formation is called the initial 
mass function (IMF). Note that / a u(M) is unequal to the gen- 
erating mass distribution /m(M), except in the random pairing 
case (see § 16. U . 

For several pairing functions we generate companions by 
drawing the mass ratio for a binary from a (generating) mass 
ratio distribution f q (q) Depending on the additional constraints 
specific for each pairing function, the resulting overall mass 
ratio distribution f q ^u{q) may or may not be equal to f q (q). The 
resulting mass ratio distribution may or may not be a function 
of primary spectral type. We refer to the specific mass ratio 
distribution for all binaries with a primary of spectral type A 
or B, for example, with f q ,ABiq)- Throughout this paper we will 
mostly use the expression (q) for the specific mass ratio 
distribution for an ensemble of binaries with a limited primary 
mass range. 

A common expression to quantify the multiplicity of a stel- 
lar population is the multiplicity fraction S (which is often re- 
ferred to as the binary fraction), defined as 



B + T + . 



S +B + T + . 



(1) 



where S is the number of single stars, B the number of binaries, 
and T the number of triple systems (e.g jReipurfh & Zinneckei 
19931) . Throughout this paper we consider only single and bi- 
nary stars, and do not consider higher-order systems, so that 
S — B/(S + B). The number of systems is given by N - S + B 
and the total number of (individual) stars is S + 2B = N(S + 1 ). 

Each pairing algorithm is provided with a (generating) bi- 
nary fraction S, which describes the fraction of stars that is 
assigned a preliminary companion star. For most pairing func- 
tions this preliminary companion is accepted as such, so that 
the overall binary fraction S d \\ equal to the generating value S. 
For several pairing functions, however, additional constraints 
are set to the properties of the companion. For example, when 
a pairing function generates a Jupiter-mass companion around 
a solar-mass star, the "primary" is usually considered to be a 
"single" star. For such pairing functions, the resulting overall 
binary fraction S a n is smaller than S. We denote the specific 
binary fraction for the set of all single stars of spectral type A/B 
and binary systems with primary spectral type A/B as Sab- 

During our analysis we make several assumptions for rea- 
sons of clarity; our models are simplifications of reality. Our 
results are not limited by these assumptions, and the models 
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Table 1. The distribution of stars over main-sequence spec- 
tral type, for the Kroupa and Salpeter mass distributions. The 
spectral type and corresponding (approximate) mass range are 
listed in columns 1 and 2 (BD = brown dwarfs). For both 
the Kroupa and Salpeter mass distribution we list the fraction 
Tn of objects and the fraction Tm of the total mass in each 
mass range. We have only considered stars in the mass range 
0.02 M < M < 20 Ma. 



SpT Mass range 



Kroupa 



Salpeter 





(M G ) 


r N (%) 


T M (%) 


T N (%) 


T M (%) 


B 


3.0-20 


1.45 


24.76 


0.10 


8.66 


A 


1.5-3.0 


2.28 


13.05 


0.18 


5.31 


F 


1.0-1.5 


2.72 


9.15 


0.22 


3.86 


G 


0.8-1.0 


2.22 


5.52 


0.18 


2.30 


K 


0.5-0.8 


7.45 


12.92 


0.61 


5.45 


M 


0.08-0.5 


51.48 


30.28 


14.09 


32.07 


BD 


0.02-0.08 


32.41 


4.33 


84.61 


42.34 


All 


0.02-20 


100.00 


100.00 


100.00 


100.00 




0.01 



0.10 



1.00 



10.00 



Mass (M s[ 



can easily be extended. Our main goal is to illustrate the impli- 
cations of adopting a particular pairing function. 

We assume that no triples or higher-order systems are 
present. Although observations have shown that a signifi- 
cant fractio n (^ 15%) of the stars a r e part of a multi- 
ple system ( Tokovinin & Smekhovl 20021; ICorreia et ai1l2 006; 
Tokovinin et al.l 20061; Eggleton & Tokovininl 120081; iHu et al.l 
2008), the properties of higher-order systems are not well un- 
derstood. Observational selection effects complicate the deriva- 
tion of the properties of these systems significantly. Higher- 
order systems are often ignored in N-body simulations due t o 



Fig. 1. The cumulative mass distributions (top) and mass dis- 
tributions (bottom) derived by Kroupa (Eq. |2j solid curve) and 
Salpeter (Eq. [3} dashed curve). The masses at which the slope 
of the Kroupa mass distribution changes are indicated with the 
vertical dotted lines. 



Chabrier 2003). In our analysis we consider two mass distri- 
butions: the Kroupa mass distribution and the Salpeter mass 
distribution. The main difference between these is the presence 
or absence of a turn over in the low -mass regime. The mass dis- 
tribution derived by lKrou pa (2001) is given by 



computational complications (e.g.. Ivan den Berk et al. 2007 ). /Kroupa(^) cc 



A full understanding of star formation, however, ultimately 
requires a full knowledge of the formation and evolution of 
higher-order systems. 

In our models the generating binary fraction for the popu- 
lation can be described with a single value S. For most pairing 
functions this results in a specific binary fraction Sm, (Mi) that 
is independent of primary mass Mi . However, in §|6] we de- 
scribe several cases where, as a result of the pairing properties, 
the binary fraction is a function of primary mass, even though 
this dependency is not included explicitly. 

Selection effects play a major role in the interpretation of 
binary star observa tions. A detailed description of the selection 
effects, such as in lKouwenhoven et alj d2007bl) . is necessary to 
derive the pairing function, the mass ratio distribution and the 
binary fraction from observations. A major bias is generally in- 
troduced by studying the binary population in a certain primary 
mass range; we describe this effect in detail for the different 
pairing functions. Throughout most of this paper we ignore the 
other selection effects. 

3. The mass distribution 

The mass distribution /m(M) defines the spectrum of masses 
in a stellar population, and is usually e xpressed as a single- 
component po wer law ( Salpeter! 1 1955b . a multi-component 
power law (e.g., lKroupall200ll) or a Gaussian distribution (e.g , 




0.01 <M/M <0.08 
0.08 < Mj M < 0.5 
0.5 <M/M Q <oo 



(2) 



The classical Salpeter mass distribution (Salpeter 119551) is 
given by 



/salpeter(M) OC M 



-2.35 



M>lM f 



(3) 



The Salpeter mass distribution is derived for intermediate-mass 
stars in the Galactic field. Although it is known that the Salpeter 
mass distribution is incorrect for masses below ~ 1 M , we 
use this mass distribution for comparison, to illustrate the effect 
of the slope of the mass distribution for low-mass stars. Note 
that, although, Eq. (f2]i introduces a turnover in Fig. Q] which 
is logarithmically plotted, there is no real turnover in the mass 
distribution (see also Appendix lA.1.3l >. 

Table Q] lists the fraction of stars of a given spectral type 
for both mass distributions. The corresponding probability dis- 

rM 

tributions /m(M) and cumulative distributions fM(M')dM 
are shown in Fig. [T| Most young stellar populations in our 
Galaxy are described accur ately with the Kroupa mass distri- 
bution dBonnell et al.ll2007h . 

In our simulations we draw objects from the mass distri- 
bution in the range 0.02 M < M < 20 M , i.e., stellar 
and substellar-mass objects. Objects less massive than M m ; n 
0.02 M are considered to be planets, and form in a different 



way than stars or brown dw arfs (see, e.g.. lPollack et al.lll996 



Kouwenhoven et al 



mass is of order M maXja b s 



2007a). The absolut e maximum stellar 
150 M Q dZinnecker & Yorke 
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20071 and references therein). The most massive star in a clus- 
ter may well depend on the total mass of the star cluster 
( Wei dner & Kroupa 2006). Stars more massive than 20 M , 
however, are extremely rare, are very short-lived, and possi- 
bly e ven form by a different m echanism than most other stars 
(e.g., Zinnecker & Yorke 2007). 



The mass distribution of 



Kroupal (1200 ll) is for all stars in 



a population, including single stars, primaries and compan- 
ions. We note that choosing primaries from an IMF and then 
choosing secondaries from a mass ratio distribution (§ 16.31 ) will 
not recover the original IMF. Thus the primary mass distribu- 
tion functi on cannot be exactly the same as the desired IMF 
(see § r7TT1 iMar Kov & ZinneckeihooU iGoodwin et al.ll2008b . 
Technically, the mass distribution should therefore not be used 
to generate, for example, a primary mass distribution. For sim- 
plicity, however, we adopt the Kroupa mass distribution as the 
generating mass distribution for each pairing function. Ideally, 
one should iteratively determine the generating mass distribu- 
tion by comparing the outcome of the pairing process with the 
Kroupa mass distribution. As this is computationally a very ex- 
pensive exercise, we skip the iteration, and simply adopt the 
Kroupa mass distribution as the generating mass distribution. 
For a proper analysis of real observations, one should keep this 
issue in mind. 



4. The mass ratio distribution 

The mass ratio distribution describes the distribution over mass 
ratio q = M2/M1 for a population of binary systems. The 
mass ratio distribution for binary systems has been s tudied 



thoroughly over t he last decades (see, e.g., Zinnecker 1984; 
Hogevee ril ll992al: T 



Mazeh et al. 2003; Halbwachs et al 



2003) 



In this paper we discriminate between three different types of 
mass ratio distributions. The generating mass ratio distribution 
is an input distribution that is used by most pairing function 
algorithms to generate binaries, although some pairing func- 
tions (e.g., random pairing) do not require a generating mass 
ratio distribution. Note that the mass ratio distribution is not 
the same as the pairing function; the mass ratio distribution is 
a property of several (not all) pairing functions; see §|6]for de- 
tails. The overall mass ratio distribution is the mass ratio distri- 
bution resulting from the pairing mechanism, for all binaries in 
the population. The specific mass ratio distribution is that for a 
sample of stars with primaries in a given mass range. The lat- 
ter is measured in observations, as a binarity survey is in most 
cases focused on a particular set of targets with given spectral 
types. 

The (specific) mass ratio distribution is usually obtained by 
a fit to the observed mass ratio distribution of the sample. The 
observed distribution is often described with a simple power- 
law: 



fr a (q) x q y " for qo < q < i 



(4) 



where the exponent y„ is fit ted (e.g. Hogeveenl 1992b; 



Shatskv & Tokovininl 120021: iKobulnicky & Fryer! 12007 



Kouwenhoven et al.ll2007b ). Distributions with y q = are flat 



while those with y q < and y q > are falling and rising with 
increasing q, respectively. Usually the adopted minimum value 



for the fit qo is the value of q below which the observations 
become incomplete due to selection effects. For distributions 
with y q < -1, a minimum value qo > is necessary such that 
f q (q) can be normalised. Sometimes the necessity of qo > is 
avoided by fitting a mass ratio distribution of the form 

/ r,(g)oc(l+g) r ' for 0<g<l (5) 

to the data (e.g. iKuiperl 19351 van Rensbergen et al .120061) . This 
distribution is more commonly used to describe the mass ratio 
distribution of high-mass spectroscopic binaries, while the dis- 
tribution in Eq. (|4]i is often used for visual binaries. As Eq. © 
and (0 show a similar behaviour (both are either falling, flat, 
or rising), we will only consider the distribution of Eq. (@). 

Alternatively, to allow for a mass ratio distribution with a 
peak in the range < q < 1, we also consider the Gaussian 
mass ratio distribution: 

(q ~ Vq) 2 ' 



/caus S (<7) oc exp < - 



2cr„ 



for < q < 1 , (6) 



where p. q and cr q are free parameters, corresponding to the 
mean and standard deviation of a Gaussian distribution with- 
out the imposed limits on q. Models with p q < and p q > 1 
show a distribution /Gauss(g) with an exclusively negative and 
positive slope, respectively. In the case where fx q «: and 
fiq » 1, the distribution may be approximated by a power- 
law (Eq.|4]l. Models with cr q » 0.5 can be approximated with 
a flat mass ratio distribution. Note that the values of p q and cr q 
do not necessarily reflect physical properties. A value p q <s 0, 
for example, merely means that the mass ratio distribution in 
the interval < q < 1 can be described by Eq. (O in this 
interval. A Gaussian mass ra tio distributions was reported by, 
Duq uennoy & Mayor! dl991 ) for solar-type stars in the solar 
neighbourh ood, who find, based on the work on the initial mass 



function by Kroup a et al. ( 1990l) . p q = 0.23 and cr q = 0.42 



For our default model we adopt a (generating) mass ratio 
distribution of the form f q (q) = 1 with < q < 1, the flat mass 
ratio distribution. 

5. The origin of the pairing function 

We define the pairing function as the algorithm which is used to 
pair individual stars into binary systems. A well-known pairing 
function is random pairing from the mass distribution. Others 
include, for example, a fixed mass ratio distribution. The pair- 
ing of binary components in a stellar population results from 
the combined effect of star formation, stellar evolution, binary 
evolution, and dynamical interactions. By studying the pairing 
of binary stars, the contributions of the latter three effects can 
be evaluated, and an estimate for the primordial binary popu- 
lation can be obtained. This primordial pairing function allows 
us to constrain the process of star formation. 

It is worth considering what pairing we might expect from 
the actual star formation process, as opposed to the various the- 
oretical constructs we describe in this paper. 

Often random pairing has been used to construct binary 
systems for various models. Random pairing has the obvi- 
ous advantage that the chosen IMF is, by design, automat- 
ically recovered. However, there is no good theoretical rea- 
son to suppose that the star formation process would produce 
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a randomly paired distribution. Furthermore, random pairing 
is ruled out observationally. The obser ved mass ratio distri- 
bution among intermediate mas s stars dShatskv & Tokovinin 
l2002t Kouwenho ven et al.l 120051 l2007allbl) and brown dwarfs 
dKraus etalJ r2008) in the nearby OB association Scorpius- 
Centaurus have indicated that the binary components are not 
randomly paired from the mass distributio n. The same result 
is fou nd for the Cygnus OB2 association dKobulnicky & Fryen 
20071) . and also by Weidner, Kroupa & Maschberger (in prep). 



will be more able to accrete material (see Goodwin et a 
and references therein; also Del gado Donate & Clarke! 



112007 



2005). 



Random pairing is further excluded by the large prevalence 
of massive binaries with a mass ratio close to unity, often 
referred to as the "twin peak" in the mass ratio distribution 
dLucy & Riccol Il979i iTokovininl l2000t iGarcia & Mermilliodl 
2001; Pinsonneaul t & Stanekl 120061: iLucvl 120061; ISoderhielml 
20071 see also §1911 

Simulations have shown that it is impossible to form sig- 
nificant numbers of binary syste ms from an initially single 
star distribution dKroupal ll995a). therefore stars in binary 
systems must have predominantly formed in binary systems. 
Observations of pre-main sequence stars also suggest that they 
have a higher multiplicity than field stars (at least for > 1 M Q ), 
suggesting th at most stars form in bina ry (or higher-order) 
systems (e.g., Goodwin & Kroupa 20051) . It is supposed that 
a primordial population with a multiplicity of almost 100% 
evolves into a field-like binary population through (a) the de- 
cay o f higher-order multiple systems (e.g. jGoodwin & Kroupa 
120051) . and (b) the dynamical destr uction of binari es in binary- 
binary encounters in clusters (e.g., Kroupafl995allcl) . Thus, the 
currently observed binary population is a complex mixture of 
primordial binaries (i.e. in the same dynamical state as when 
they formed), and dynamically evolved binaries (which may 
have different characteristics, or even companions to their ini- 
tial state). 

Simulations of binary star formation have comprehensively 
failed to produce systems that match observations, even when 
the dynamical evolution o f the initial states is accounted for 
(see iGoodwin et al.l 120071 and references therein). However, 
hydrodynamic simulations of star formation suggest that com- 
panions usually form by the fragmentation of massive, disc- 
like circumstellar a ccretion regions around young stars (see 



Goodwin et al.ll2007 and references therein). 



In such a scenario for companion formation it would be ex- 
pected that the secondary should have a roughly similar (i.e. 
within a factor of three or four) mass to the primary, espe- 
cially at low separations. A massive enough region to frag- 
ment is only present during the earliest (e.g. class 0/1) phases 
of star formation before the star(s) have accreted the majority 
of their natal core. Thus, the secondary will be present whilst 
a large reservoir of gas is also present around it. In the case 
of a star that will eventually grow to be (for example) 5 M 
the secondary will form whilst the primary is only 0(1 M ) 
and several solar masses of gas are present in the circum- 
stellar environment. A secondary will presumably form with 
an initial mass close to the opacity limit for fragmentation, 
~ 10~ 2 M . However, it is difficult to imagine a scenario in 
which the secondary fails to accrete at least some of the circum- 
stellar material, especially as the secondary will form with an- 
gular momentum similar to that of the accretion region, and so 



Therefore, the secondary mass is expected to be a reasonably 
large fraction of the primary mass. In particular, it should be 
difficult for a companion to a B-star to remain at brown dwarf 
or M-star mass due to the large amount of material available for 
accretion. In particular, we would expect a rough correlation 
between separation and mass ratio, with closer co mpanions be 



ing g enerally more massive as is observed (e.g. iMazeh et al 
Il992l for field G-dwarfs). 

In addition, dynamical evolution in clusters will act to 
destroy the most weakly bound systems (i.e. the widest and 
lowest-mass companions), further biasing the mass ratio distri- 
butions away from low-g. 

We would therefore argue that random pairing over the full 
mass range is the last type of pairing that would be expected 
from the star formation process (see, however, § 16.1 .It . 

6. Analysis of frequently used pairing functions 

There are many ways to obtain a population of binary sys- 
tems from a mass distribution /m(M). We analyse the most fre- 
quently used algorithms in the sections below. In general, the 
masses of the members of a stellar grouping are drawn from the 
mass distribution. A fraction of the stars is assigned a compan- 
ion star (either from the mass distribution, or using a mass ratio 
distribution), or the mass is split into a primary and companion 
star. Four commonly used mechanisms are the following: 

- Random pairing (RP). The masses of both binary com- 
ponents are randomly drawn from the mass distribution 
/m(M). For each system, the most massive component is 
labelled "primary star", the other component "companion 
star" (see § 16. 11 1. 

- Primary-constrained random pairing (PCRP). The primary 
mass M[ is drawn from the mass distribution /m(M), and 
the companion mass M2 is chosen from the same mass dis- 
tribution, with the additional constraint that M2 < Mi (see 



- Primary-constrained pairing (PCP). The primary mass is 
drawn from the mass distribution /m(M). The companion 
mass is then determined by a mass ratio that is drawn from 
a distribution f q {q), with < q < 1 (see below). 

- Split-core pairing (SCP). The total mass of the binary is 
drawn from the mass distribution /m(M). The mass ratio of 
the binary is then determined by a mass ratio distribution 
which is drawn from f q (q) with < q < 1. Finally, the 
primary mass M; = M(l + q)' 1 and companion mass M2 = 
M(l + q~ l )~ l are determined (see below). 

In the case of PCP and SCP there is another complication, 
which occurs if M2 = qM\ < M2, m m- In this case, the mass ratio 
distribution generates a companion mass smaller than the per- 
mitted value M2, m i n , for example the deuterium-burning limit, 
while such objects are usually not considered as companions 
in a binary system. There are three straightforward choices on 
how to handle such companions: 

- Accept all companions (PCP-I/SCP-I). All companions are 
accepted, regardless of their mass. The resulting mass ratio 
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distribution obtained with this method is equal to the gen- 
erating mass ratio distribution (see § 16.3. H and § 16.4. 11 1. 

- Reject low-mass companions (PCP-II/SCP-II). If 
^2,min/.Mi,max < q m i n , a fraction of the companions 
has M2 < M2 y m m - All companions with M2 < M2, m i n are 
rejected, and the corresponding "primaries" are classified 
as single stars (see § 16.3.21 § 16.4.21 and below). 

- Redraw low-mass companions (PCP-III/SCP-III). For all 
binaries with a companion mass M2 < A^min the mass ra- 
tio is redrawn from f q {q). This procedure is repeated until 
M2 > M2 im in. This method is equivalent to drawing a mass 
ratio from the distribution f q {q) with limits M2,m m IM\ < 
q < 1 (see § IQ31 § l6A3l and below). 

It is possible in pairing function SCP-II that the resulting single 
star mass is smaller than M m ; n , and for SCP-III the splitting up 
is not possible if the binary system mass is smaller than 2M m ; n . 
In these cases there are three possibilities of dealing with this 
problem: 

- Accept all singles/primaries (SCP-IIa/SCP-IIIa). If a pair- 
ing mechanism produces a single or primary star less mas- 
sive than M2, m in, it is accepted, and included in the model. 

- Reject low-mass singles/primaries (SCP-IIb/SCP-IIIb). If a 
pairing mechanism produces a single or primary star less 
massive than Mijmn, it is rejected and removed from the 
model. 

- Do not split-up low-mass cores (SCP-IIc/SCP-IIIc). Cores 
with a mass Mj < 2M m ; n are not split up; these become 
single stars of mass M$ - Mj. 

For reasons of simplicity and clarity, we use a generating 
mass distribution with a minimum value of 2M m ; n for the SCP 
pairing functions, and adopt a minimum mass M m j n for the 
companions. Implicitly, we therefore only consider the vari- 
ants SCP-IIb and SCP-IIIb in this paper. From hereon, we use 
"SCP-II" and "SCP-III" to refer to SCP-IIb and SCP-IIIb, re- 
spectively. 

The above-mentioned pairing functions (RP, PCRP, and the 
three variations of PCP and SCP) are described in detail in the 
subsections below, while their differences are discussed in §|7] 
Unless stated otherwise, we have adopted a Kroupa generat- 
ing mass distribution, a flat generating mass ratio distribution 
(where applicable), and a generating binary fraction S of 100% 
(see §[8] for a discussion of these assumptions). The main dif- 
ferences between the (renormalised) resulting mass ratio dis- 
tributions are show in Fig. [2] where the top panel represents 
the overall mass ratio distribution. The middle and bottom pan- 
els represent the specific mass ratio distributions f q \M x (q) (the 
subscript M\ indicates a restricted primary mass range) for bi- 
naries with high-mass primaries and low-mass primaries, re- 
spectively. Note that the derived mass ratio distribution for a 
sample of stars does not only depend strongly on the pairing 
function, but also on the targeted sample of stars. For the same 
reason, the (renormalised) companion mass distribution, shown 
in Fig. [3] depends strongly on the primary mass range. A two- 
dimensional version of Fig. [2]is shown in Fig. [4] In the sections 
below we discuss in detail the pairing functions and the above- 
mentioned figures. 



Note that the choices made in this paper do not imply 
that stellar populations indeed have these properties. It is not 
known how binary stars are formed, so that no robust pre- 
dictions of their properties can be made. Different binary for- 
mation mechanisms may produce different mass ratio distri- 
butions, possib ly varying with primary mass, period, or ec- 
centricity (e.g. | Heggid Il975c iKrumholz & Thompson! [2007; 
IZinnecker & Yorka 120071) . In addition, dynamical evolution 
after the formation process may alter the bin ary fraction 



Hills 1975 



and the mass and mass ratio distributions (e.g. 
i Heggidl 19751) possibly as a function of environment (see, e.g. 
iKroupa et alH"l999: Preibisch et alJl2003l buchene et al]l"2004 



iKohler et alJl2006l:lReipurth et alJl2007L and §|92]l. Other pair- 
ing functions suggested in li terature includ e ra ndom pairing 
over a restricted mass range dKroupall 1 995aH3lbl: IKroupa et alj 
12003b iThies & Kroupijl2007i see also § |6XTI gravitationally- 
focused random pairing (Kouwenhoven et al, in prep.), ordered 
pairing (Oh et al, in prep. ), pairing resul ting from the disso- 
lution of small-jV clu sters dClarkelll996bl) . "two-step" pairing 
( Durisen et alj 200 binary for mation from ring fragmenta- 
tion (iHubber & Whitworfh||2005b . and numerous others. 



6.1. Pairing function RP (random pairing) 

In the case of random pairing (RP), the primary and com- 
panion mass are both independently drawn from /m(M), and 
swapped, if necessary, so that the most massive star is the pri- 
mary. As a result of this swapping, neither the resulting pri- 
mary mass distribution ,/m,(Mi), nor the companion mass dis- 
tribution /^(Mj), nor the system mass distribution /m t (Mt) 
is equal to the g enera t ing ma s s distribution ./m(M); s e e, e,g ., 
Warner! dl96lh : iToutl d 199 lb ; iMal'Kov & Zinneckerj d200lb . 
On the other hand, the mass distribution of all stars / a ii(M), 
i.e., all singles, primar i es and companions, is equal to fm(M). 
MaFKov & Zinneckerl (120011) derived general expressions for 
the distribution over primary star mass Mi, companion star 
mass M2, and system mass Mj — M\ + M2, respectively: 

f Ml (M 1 )= 2/ M (Mi) j c Ml f M (M) dM (7) 

/m 2 (M 2 )= 2f u {M 2 )faj M (M)dM (8) 

fu T (Mr) = 2 J c Mt ~ c f M {M) f M (M T - M) dM (9) 

Below we calculate the mass ratio distributions resulting 
from random pairing. Let the generating mass distribution / be 
defined on the interval [c, d]. For random pairing we draw two 
stars with masses x and y from /, and swap, if necessary, so 
that the primary is the most massive star. The resulting mass 
ratio distribution is bounded by < c/d < q < 1. To derive 
the overall mass ratio distrib ution we follow the appendix of 
Piskunov & MaFKovl dl99lb and derive f q (q) from its cumula- 
tive distribution function F q (q) = P(y/x < q): 

F q (q)= ff U{x,y)dxdy = f C 2f(x)f(y)dydx, (10) 

where the factor 2 accounts for the fact that the pairs of masses 
are swapped in order to ensure that y < x. The integration do- 
main is shown in the left-hand panel of Fig. [5] The probability 
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Fig. 2. The (renormalised) mass ratio distributions resulting from the different pairing functions. From top to bottom, the panels 
show the overall mass ratio distribution (i.e., for all binaries), the mass ratio distribution for binaries with 1 .5 M Q < M\ < 20 M , 
and for binaries with 0.02 M < M\ < 0.08 M . The models consist ofN = S+ B = 5x 10 5 particles and all have a generating 
binary fraction ofB- 100%. For each model we adopt a Kroupa generating mass distribution in the mass range 0.02 - 20 M Q , 
and, when applicable, a flat generating mass ratio distribution f q (q) = 1 (0 < q < 1). This figure illustrates that each pairing 
function results in a different overall or specific mass ratio distribution (see also Figs.[8]and|7]i. 
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Fig. 3. The (renormalised) companion mass distributions resulting from the different pairing functions, for the models shown in 
Fig. [2] The curves indicate the distribution for all binaries (dashed curves), for binaries with 1.5 M < Mi < 20 M (dotted 
curves) and for binaries with 0.02 M < M\ < 0.08 M (solid curves). 



density f q for q (i.e., the overall mass ratio distribution) is then 
given by the derivative of Eq. (fTOb . 

For observational reasons, surveys for binarity are often re- 
stricted to a certain range of primary spectral types. To derive 
the effects of the selection on primary mass, the derivation of 
fq(q) again proceeds via the cumulative distribution function, 
which is now given by: 



F q (q) = P(y/x < q\x\ < x < xi) oc JJ f lp (x,y)dxdy , 

S' 



(11) 



where the primary mass range is restricted to the range [x\, x^\. 
The integration domain 5" is as shown in Fig. The inte- 
gration limits now depend on whether q is larger or smaller 
than c/x\. For q < c/x\ the integration domain S' is given by 
c/q<x<X2/\c<y<qx (middle panel in Fig. [5j, while 
for q > cjx\ S' is defined by x\ <x<X2Ac<y<qx 
(right panel in Fig. 01. The minimum possible value of q is 
c/x2, so that F{q) = for q < c/x2- The general expression for 
the cumulative mass ratio distribution for a sample of binaries 
with a restricted primary mass range is given in Appendix lA. 1 1 
The dependence of the mass ratio distribution on the primary 
mass range for RP is shown in Fig. |6]for the Kroupa (top) and 



Salpeter (bottom) generating mass distributions. Note that each 
mass ratio distribution is renormalised such that its maximum 
is unity. Note the differences between the overall mass ratio dis- 
tributions fq.aiiiq) in the left-hand panels of Fig. [6] The Salpeter 
mass distribution results in a peak in the overall mass ratio dis- 
tribution at q ss 1 . The Kroupa mass distribution, on the other 
hand, is on average much shallower, and therefore produces a 
peak at small q (cf. the solid lines in the top panels of Fig. lATl i. 

In the special case where the mass distribution is of the 
form /m{M) oc M~ a with a + 1, the overall mass ratio dis- 
tribution is given by 

fM) - y 



(c/of) 2r 



q 



(12) 



(i - {cidy? 

where y = a — 1. Realistic mass distributions cover a broad 
range of masses, i.e., c <K d, for which the expression simplifies 
to: 

'(a - l)q a ~ 2 a > 1 
(1 - a)q~ a a < 1 

(Pis kunov & MarKovlll99lh . For distributions with a w 1, the 
above expression is not a good approximation, and Eq. $Y2\ 
should be used. In the case of a power-law mass ratio distribu- 
tion with c <s d, the overall mass ratio distribution resulting 



(13) 
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Fig. 4. The distribution of binary systems over primary mass M\ and mass ratio q for the different pairing functions. All popula- 
tions have a Kroupa generating mass distribution, a generating mass ratio distribution of the form f q (q) - 1 (when applicable), 
and a generating binary fraction S of 100%. From left to right, the panels show the results for pairing functions PCP-I, PCP- 
II, PCP-III and RP in the top row, and SCP-I, SCP-II, SCP-III and PCRP on the bottom row. From light to dark, the contours 
bound values of 0.01%, 0.02%, 0.05% (thick curve), 0.1% and 0.2% for (q, log Mi). Each model contains 10 6 binaries. Each dot 
represents a bin in mass ratio and logarithmic primary mass. No selection effects have been applied. 
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Fig. 5. The integration domains for pairing functions RP and PCRP. Left: the integration domain in the (x, y) plane (where x and y 
are the primary and companion mass, respectively) for the determination of the cumulative distribution of q — y/x. The lower and 
upper limits on the generating mass distribution are given by c and d, which implies that c/d < q < 1. The integration domain S 
is given by: c/q < x < d A c < y < qx. Middle and right: same, for the case that the primary mass x is restricted by x\ < x < Xi. 
If q < c/xi the integration domain S' is given by c/q <x<X2Ac<y<qx (middle panel), while for q > c/xi, the domain 5" is 
defined by xi <x<X2Ac<y<qx (right panel). 



from random pairing can thus be described by a simple power- 
law. Note that Eq. (TT31 does not depend on c or d; see also 
the bottom-left panel in Fig. [6] For a Salpeter mass distribu- 
tion (a = 2.35), for example, the resulting overall mass ratio 
distribution is fq A \\(q) q 035 . The general expression for /„ 
when the primary mass range is restricted to the range [xi , xj\ 
is given in Appendix lA. 1 1 For a sample of binaries with a small 
primary mass range (x\ ~ X2), fq(q) is proportional to q~ a for 



c/x\ < q < 1 and zero otherwise: 
(0 



fqMM) 



79 



q < cjx\ 
cjx\ < q < 1 



(14) 



U-(c/-»'i) T 

Random pairing from a mass distribution that is approximately 
a single power-law (e.g. the Kroupa IMF) thus generally re- 
sults in a three-segment mass ratio distribution, which is zero 
for q % c/x2 and exhibits a peak at q m c/x\ (Eq. IA. 15T >. As 
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Fig. 6. How does the specific mass ratio distribution / 9; m, (q) resulting from random pairing (RP) depend on the mass distribution 
and the selected primary mass range? This figure shows the (renormalised) resulting mass ratio distributions for the Kroupa mass 
distribution (top panels; Monte Carlo simulations) and the Salpeter mass distribution (bottom panels; analytical calculations) in 
the mass range c < Mi < 20 M . The lower limit of /m{M) is set to c = 0.02 M for the solid curves, and c = 0.08 M for the 
dashed curves; the difference between these models is thus the presence of brown dwarfs. In all panels the models with brown 
dwarfs have on average a lower mass ratio. The left-hand panel shows the overall mass ratio distribution (for all binaries in the 
population). The other panels show the specific mass ratio distribution for the binaries with a primary of the spectral type shown 
in the top-right corner of each panel (BD = systems with brown dwarf primaries). All distributions are normalized so that their 
maximum is unity. Each model contains 5 x 10 5 binaries and has a binary fraction of 100%. 
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Fig. 7. The specific binary fraction Sm,(Mi) as a function of primary mass for pairing functions RP, PCP-II, SCP-I, SCP-II and 
SCP-III. Models with a generating binary fraction S of (from bottom to top) 5%, 50% and 90% are shown. The horizontal dotted 
lines in each panel represent the generating binary fraction. We adopt a Kroupa generating mass distribution with c = 0.02 M 
and d — 20 M . For pairing functions PCP-II, SCP-I, SCP-II and SCP-III we show results for models with generating mass ratio 
distributions f q {q) oc q 7 , with y — (solid curves), y = -0.3 (dash-dotted curves), and y = -0.6 (dashed curves). 



already noted by IZinnecken d 19841) . a sample of binaries with 
high-mass primaries thus shows a peak in f q {q) at low q, while 
a sample with low-mass primaries peaks at high q. 

For pairing function RP the overall binary fraction T> d \\ is 
equal to the generating binary fraction < B. The specific binary 
fraction Sm^Mi), however, depends on the surveyed primary 
mass range: the larger the primary mass, the higher the spe- 
cific binary fraction (unless & = 100%). This can be under- 
stood as follows. For random pairing, the mass of the primary 
stars is drawn from the mass distribution. A fraction & of the 
primaries is assigned a companion, and primary and compan- 
ion are swapped, if necessary, so that the primary is the most 
massive star. This swapping leads to an increased number of 



binaries with a high-mass primary, and a decreased number of 
binaries with a low-mass primary, and hence a mass-dependent 
binary fraction. The relation between specific binary fraction 
and primary mass for random pairing is given by 



(15) 



where Fm{M\ ) is the cumulative mass distribution evaluated at 
mass Mi (see Appendix |A.1.2| i. Fig.Qshows the binary frac- 
tion as a function of primary mass for several binary fractions, 
for the Salpeter mass distribution. Clearly, we have &m, (Mi) = 
1 when & = 100%, S Ml (Mi) = when & = 0%, and 
S M] (M min ) = 0. Systems with Mi > <M> have B Ul (M{) > S, 
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while systems with M x < (M) have B Ml (M{) < S (see 
Appendix |A. 1 .21 for a full derivation). The magnitude of this 
difference depends on the shape of the mass distribution and 
the generating binary fraction. Note that the variation of binary 
fraction with primary mass is purely a result of the choice of 
pairing function; no explicit variation of binary fraction with 
primary mass is included in the simulations. 

6.1.1. Restricted random pairing (RRP) 

Restricted random pairing (RRP) is very similar to random 
pairing (RP) as described in § I6.1I with the difference that the 
binary components are now drawn from a limited mass range. 
All properties derived in § I6. H are thus applicable to the result- 
ing binary sub-population resulting from RRP. However, the 
nature of RRP implies the presence of one or more other sub- 
populations that have a formed via another process. The other 
sub-populations could have alternative pairing function, such 
as RRP with different lower and upper mass limits, or a com- 
pletely different pairing function. 

Kroupal dl995allc1lbh finds that observations of binary sys- 



range \x \ , X2] now becomes: 



tems are consistent with the population being born with pair- 
ing function RRP in the stellar mass range, prior to the effects 
of pre-main sequence eigenevolution. Further motivated by the 
difference between the observed mass ratio distribution and 
semi-major axis distribution of binary systems w ith a stellar 
primary and those with a brown dwarf binary (e.g . . iBouv et al 



2003 



Burga sser et alJ 120031: iMartfn et alJ 12003b IClose et al 



2003), this implies that t he brown dwarf popu l ation has formed 
with a differ ent process Kro upa et alJ (12003b : iThies & Kroupa 
(20071 120081) . Their proposed model with stellar and substel- 
lar sub-populations is further supported by the existence of the 
brown dwarf desert among solar-type stars (see § 19. 11 1. 

6.2. Pairing function PCRP (primary-constrained 
random pairing) 



For primary-constrained random pairing (PCRP), each primary 
mass Mi is drawn from the mass distribution /m(M) with limits 
c < M < d. The companion mass Mo is also drawn from the 
same mass distribution, but with the additional constraint that 
Mo < Mi . The limits on the resulting mass ratio distribution are 
equivalent to those of random pairing < c/d < q < 1. Writing 
the primary mass distribution as f(x) and the re-normalised 
companion mass distribution as f'(y), the expression for the 
joint probability distribution is f VCIV (x, y) = f x (x)f' y (y), which 
is normalised to unity due to the re-normalisation of f y (y). To 
derive the overall mass distribution f q one can proceed as for 
the RP case (see § 16. 11 1. The integration domain S is again as 
shown in Fig. [5] and 



F(q) = P(y/x < q) 



H rqx 

/(*)/' 

«J clq Jc 



(y)dydx . 



(16) 



F q (q) ■ 



q < c/x2 

L7q fc* fWf'Wydx c/x 2 <q< c/xi 

X 2 L qx fW(yWy dx C M ^ 9 1 



(17) 



where the integration domain 5" is as in Fig. [5] The expres- 
sion for the resulting / 9 , a ii(<?) in the case of a single power- 
law generating mass distribution is given in Appendix |A.2| For 
PCRP the distribution f q ^\\{q) contains a term that diverges for 
q t 1, which can be seen in Fig. [2] (see also Eq. (lA.22b in 
Appendix IA.21 >. 



Note that the normalisation constant for f'(y) depends on x. 
The expression for the sample with a restricted primary mass 



For a sample with primary masses restricted to the range 
[xi,X2], Eq. ( fTTI i has to be worked out. We will not explic- 
itly show the results here. Most importantly, f g (q) is zero for 
q < c/x2, and exhibits a peak at q w c/xi. The distribution 
for a sample of binaries with high-mass primaries thus peaks at 
low q, and the distribution for low-mass binaries peaks at high 
q. If the primary mass range is small (xi w X2), fqiq) can be 
approximated with Eq. ( TT4l >. For a sample of stars with a very 
small primary mass range, the mass ratio distributions result- 
ing from PCRP and RP thus give the same results. Differences 
between the two pairing functions become larger for realistic 
primary mass ranges. 

For pairing function PCRP, the companion mass distribu- 
tion Jm^Mj) for a set of primaries of identical mass is equal 
to the generating mass distribution /m(M) in the mass range 
M m j n < Mo < Mi . The companion mass distribution can thus in 
principle be used to derive the properties of the generating mass 
distribution. For example, in a set of binaries with a primary 
mass of 1 M , those with mass ratio q < 0.08 have brown dwarf 
companions. If the observations are of good enough quality 
to study the mass ratio distribution below q = 0.08, and it is 
known a-priori that the pairing function is PCRP, the mass ra- 
tio distribution can be used to constrain the slope of the mass 
distribution in the brown dwarf regime. 

The pairing algorithms RP and PCRP appear similar, but 
their difference is for example seen in the primary mass dis- 
tribution. For RP there is a larger number of binary systems 
with high-mass primaries, which can be understood as fol- 
lows. Suppose a primary mass Mi of 5 M is drawn from the 
mass distribution. For PCRP, the companion mass M2 is al- 
ways smaller than the primary mass, while for RP the com- 
panion mass can take any value permitted by the mass distri- 
bution (i.e. also M2 > Mi, after which the components are 
switched). Another difference is that, unlike RP, for PCRP 
the binary fraction is independent of primary spectral type: 
Sail = 2?m, (Mi) = S. For realistic mass distributions (e.g., 
Salpeter- or Kroupa-like), the overall mass ratio distribution 
f qA n{q) of PCRP is peaked towards high values of q, while that 
of RP is peaked towards low values of q. The pairing functions 
RP and PCRP can be excluded if more than a 1 - 2% of the 
intermediate-mass stars are 'twins' (q > 0.8). 
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6.3. Pairing function PCP (primary-constrained 
pairing) 

In models with primary-constrained pairing (PCP), each bi- 
nary system is generated by drawing a primary mass M\ from 
/m(M) in the range c < M < d, and a mass ratio q from the 
generating distribution h q {q). The companion mass is then cal- 
culated from M2 = qM\ . Due to the nature of this pairing mech- 
anism it is possible that the resulting companion is of very low 
mass, for example a planetary mass if a very small mass ratio 
is drawn. Below we describe three variants of pairing function 
PCP, each of which handles very low mass companions in a 
different way: accepting all companions (PCP-I), rejecting the 
very low-mass companions (PCP-II), and redrawing the mass 
ratio if the companion mass is of very low mass (PCP-III). 

6.3.1. Pairing function PCP-I 

PCP-I is the simplest variant of PCP: the primary mass is drawn 
from /m(M) and the mass ratio from h q (q), and no further con- 
straints are set. As a result, the specific mass ratio distribution 
fq-,M\ (?) an d overall mass ratio distribution f q #\\(q) are equal to 
the generating mass ratio distribution h q (q). Additionally, the 
specific binary fraction < BmSM\) and overall binary fraction 
S a n are equal to the generating binary fraction S. 

The companion mass for PCP-I can be arbitrarily small: 
A^2,min = 0. Several companions may thus have masses signif- 
icantly lower than the deuterium burning limit (~ 0.02 M Q ). 
Even planetary companions are considered as companion 
"stars" for the pairing function PCP-I. However, if we do in- 
clude planets, we make the implicit assumption that the star for- 
mation process is scalable down to planetary masses. This as- 
sumption is in contradiction with the theories t hat suggest that 



tion is given by: 



stars a nd brown dwarfs form b y fragmentation dGoodwin et al 
2007t IWhitworth et"ai] 120071). while planets form by core- 
accretion (see, e.g.. IPollack et al.lll996T) . This is an important 
point to keep in mind when using PCP-I, i.e., when adopting 
a mass ratio distribution that is fully independent of primary 
spectral type. 

6.3.2. Pairing function PCP-II 

For PCP-II, companions with M2 < c, are rejected, and the cor- 
responding primary stars are classified as single stars. There are 
two reasons why one may want to consider using PCP-II. First, 
one may wish to use this prescription if a minimum compan- 
ion mass is expected from theory, for exam ple the Jeans mass 
or the opacity limit for fragmentation (c t> IHovIpIIiqs^HRrpsI 



1976; 



1982; 



Low & Lvnden-Bell fll 

Larsonlll969l. ll992ll20"05l:lMasunaga & Inutsukall2000l) . 



The second, more observational approach may be to "ignore" 
the low-mass companions. Although this seems somewhat ar- 
tificial, this method is often used in practice. Planets are usu- 
ally not considered as companions (the Sun is a "single star"), 
which implies a limit c = 0.01 - 0.02 M . 

Due to the rejection of low-mass companions the overall 
mass ratio distribution is zero for < q < c/d. In the range 
c/d < q < 1, the expression for the overall mass ratio distribu- 



fqMlil) 



Jc/q 



f pcp (M u q)dM ] 



= k f h q 

Jc/q 



(q)f Ml (M 1 )dM l , 
(18) 

where h q (q) is the generating mass ratio distribution and k is 

a normalisation constant which ensures that J Q f q ,a\\(q)dq = 1. 
The distribution f q .s.\\{q) has a higher average mass ratio than 
the generating mass ratio distribution h q (q) as a result of re- 
jecting the low-mass companions. For a a sample with a re- 
stricted primary mass range, c < x\ < M\ < xt < d, the 
expression for the specific mass ratio distribution is given in 
Appendix lA.3l by Eq. (IA.29b . The specific mass ratio distribu- 
tion is given by a three-segment powerlaw, with slope changes 
at q — c/x2 and q = cfx\. Note that for q > c/x\, the distribu- 
tion f^M, (q) is equal to the generating mass ratio distribution. 
For a sample of high-mass primaries, where c/x\ < 1, we thus 
have (q) ~ h q (q). For a sample of low-mass primaries, 
c/x2 < q < c/x[ for most of the mass ratio range, and the cor- 
responding term dominates. For these binaries, / ?; m 1 (?) differs 
significantly from h q (q). 

The specific binary fraction for a sample of systems with 
primary mass Mi is given by 



(M0 = S 



f 



h q (q) dq < S . 



(19) 



(Mi) 



Note that B Ml (Mi) is independent of f M (M) for PCP-II. Fig. fj] 
shows the specific binary fraction Sm,(Mi) as a function of S 
and f q {q). The specific binary fraction depends on the shape of 
the mass ratio distribution, and is independent of the mass dis- 
tribution. The overall binary fraction T>. d \\ after rejection of the 
low-mass companion is smaller than the generating binary frac- 
tion £ (see Eq. (IA.321 > in Appendix IA. 3. 2l i. The binaries with 
high-mass primaries are hardly affected by the rejection algo- 
rithm. For these binaries the specific mass ratio distribution and 
specific binary fraction are practically equal to those for PCP-I: 
the specific binary fraction is equal to S, and the specific mass 
ratio distribution is equal to h q {q). For the very low-mass pri- 
maries, however, a large fraction of the companions is rejected, 
and therefore the specific binary fraction is low. The remaining 
companions of these stars have a mass comparable to that of 
their primary, and the resulting mass ratio distribution for the 
lowest-mass binaries is peaked to unity (see Fig. |2j. 

As an example, consider a stellar population with S = 
100% and a generating mass ratio distribution f q {q) = 1 and 
A^min - 0-02 M©, systems with B-type primaries and systems 
with M-type primaries have a resulting specific binary fraction 
of 99% and 87%, respectively. If we also do not consider brown 
dwarfs as companions (so M2, m i n = 0.08 M ), then the specific 
binary fractions are 96% and 51%, respectively. 

6.3.3. Pairing function PCP-III 

For PCP-III the primary mass is drawn from /m(M) in the range 
c < M\ < d, and the mass ratio is drawn from h q (q). If the re- 
sulting companion star mass is smaller than c, the mass ratio 
is redrawn from h q (q) until a companion with mass M2 > c is 
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obtained. This is equivalent to renormalising h q (q) in the range 
<7min(^i) < q < 1, where q m m(Mi) = c/M^ This effectively 
results in a mass-dependent generating mass ratio distribution 
h'q(q). The expression for the resulting overall mass ratio dis- 
tribution f qA \\{_q) are then: 

f q Mq)= I fpcp(M u q)dMi = I h' q (q)f Ml {Mi)dMi. 

Jc/q Jc/q 

(20) 

The specific mass ratio distribution for the restricted primary 
mass range [x\, is given in Appendix |A. 31 A sample of bi- 
naries with high-mass primaries has c/x\ <K 1. For high-mass 
binaries h' q (q) as h q (q), and therefore fq-M,(q) ~ h q (q). For the 
very low-mass primaries, however, a large fraction of the com- 
panions is redrawn. Consequently, all binaries with a low-mass 
primary have a mass ratio close to unity. The resulting mass 
ratio distribution for the lowest-mass binaries is thus peaked to 
unity (see Fig. [2j. 

As a result of the redrawing of the companions for pair- 
ing function PCP-III, the resulting overall binary fraction S a n 
equals the generating binary fraction S, and the specific binary 
fraction Sm, (Mi) equals & for any primary mass range. 

6.4. Pairing function SCP (split-core pairing) 

For split-core pairing (SCP) one assumes that the system "core" 
mass Mc is drawn from a core mass distribution /m c (Mc) with 
2ec < Mc < 2ed, where e is the star forming efficiency. Split- 
core pairing is frequently inferred from observations of dense 
cores in star forming regions, assuming that a fraction of the 
cores fragment into bin aries (see, e.g.. iGoodwin et al.l 12008; 
Swift & Williamsl l2008). As a core collapses, it forms one or 
two stars with a total mass Mj = eMc- The resulting minimum 
and maximum primary masses are thus c and 2d, respectively. 
The star forming efficiency may be a function of various pa- 
rameters, for example the mass of the core. For simplicity in 
our analysis, however, we keep the star forming efficiency fixed 
to e = 1 for all values of Mc- The total mass of each binary is 
thus Mt = Mc- Note that "random fragmentation" (random 
splitup of a clump into two stellar components) is very differ- 
ent from random pairing of two components from the IMF (see, 
e.g., Figs.|2]and[3]l. 

For pairing function SCP, the binary total mass is thus 
drawn from a distribution /m c (Mc)- Note that, although we 
adopt /m c {M c ) = fKmupi(M) in this paper, there is no obvi- 
ous prerequisite that /m c (Mc) should be a standard IMF. The 
binary is split up according to a mass ratio that is drawn from 
a generating mass ratio distribution h q (q). Given the core mass 
Mc — Mi + M2 and the mass ratio q, the primary and compan- 
ion mass are given by 



M c 
q+l 



and M2 



Mr 



+ 1 



(21) 



respectively. As a result of this procedure, it may happen that 
a companion mass smaller than the minimum mass c is drawn. 
Similar to pairing function PCP (§ 16.31 ), there are three ways to 
address this issue: accepting the low-mass companions (SCP- 
I), rejecting the low-mass companions (SCP-II), and redrawing 



the mass ratio if the companion mass is too low (SCP-III). We 
discuss these three variants of SCP in the sections below. 



6.4.1. Pairing function SCP-I 

For pairing function SCP-I, all binary components resulting 
from the split-up mechanism are accepted, irrespective of their 
mass. Stars with substellar and planetary companions are thus 
also considered as "binary stars". Due to the nature of this 
pairing process, the overall mass ratio distribution is equal to 
the generating mass ratio distribution. The specific mass ra- 
tio distribution, however, is a function of spectral type. A full 
derivation of the speci fic mas s ratio distribution is given in 
Appendix I A.4I see also IClarkd(ll996ah . In most cases the pri- 
mary mass range is contained within the range [2c, d], which 
corresponds to "case 7" in Appendix lA.4l Consider the special 
case of a single power-law mass distribution /m c (Mc) 00 M c a 
(a + 1) and a uniform mass ratio distribution h q (q) = 1. Under 
these assumptions, the expression for the specific mass ratio 
distribution is fi f M,(q) K (1 + q) l a for < q < 1, if either 
the primary mass range is contained within the range [2c, d] 
(case 7), or if x\ = Xz- Note that this expression is identical to 
that in Eq. ©. The highest mass binaries thus have on average a 
low mass ratio, and the lowest mass binaries have on average a 
high mass ratio. Note that these trends are present, even though 
the generating mass ratio distribution produces mass ratios in 
the range < q < 1, irrespective of the core mass Mr- 

Pairing function SCP-I naturally results in a mass- 
dependent binary fraction: 

S/m,(M,) 

&Mi(M\) - — — — — — — — — — — — — , (22) 

(Mi) + (1 - S)f Mc (Mj ) 

where /m,(Mi) is the primary mass distribution, /m c (Mi) the 
generating (core) mass distribution evaluated at mass Mi , and 
S the generating binary fraction; see Appendix IA.4.21 for the 
derivation. The mass-dependence of the binary fraction for 
SCP-I is illustrated in the middle panel of Fig. [7] In general, 
S = 100% for c < Mi < 2c, as /m c — in this mass range. 
For 2c < M; < d the specific binary fraction is more or less 
independent of Mi , while beyond Mi = d, it decreases down to 
zero at Mi = 2d. The latter dependence is due to the fact that 
the cores are more massive than the primary stars they poten- 
tially form. As a result, the high-mass targets are dominated by 
cores that have not split up. The overall binary fraction S d n is 
always equal to S. 

6.4.2. Pairing function SCP-II 

For SCP-II the companion is rejected if M2 < c, and the pri- 
mary star becomes single. The resulting primary mass range is 
then c < Mj < 2d - c, and the companion mass is in the range 
c < Mi < d. The full derivation for the overall and specific 
mass ratio distributions is given in Appendix I A.4I In most re- 
alistic cases the primary mass range [jci, JC2] is fully enclosed 
in the mass range 2c < Mi < d ("case 1"). Unlike SCP-I, the 
resulting overall mass ratio distribution is unequal to the gener- 
ating mass ratio distribution, but contains more high-g binaries 
instead. 
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The specific binary fraction B Ml (Mi) resulting from SCP-II 
varies with primary mass M\ : 



gioo(Mi)g/M,(Mi) 
Sf Ml (Mi) + (l-S)fM c (Mi) 



(23) 



where /m, is the primary mass distribution, /m c is the gen- 
erating (core) mass distribution, and & is the generating bi- 
nary fraction. The quantity Sioo(Mi) represents the value of 
S Ml (Mi) for a population with S = 100%. The full derivation 
of Eq. d23l is given in Appendix |A.4.2| 

Fig. |7] shows the specific binary fraction S(Mj) as a func- 
tion of S and h q {q). The specific binary fraction decreases with 
decreasing primary mass, as, on average, more low-mass com- 
panions are rejected among lower-mass primaries. The major- 
ity of the newly formed single stars (due to rejection of low- 
mass companions) is thus of most low mass. For the very 
lowest-mass stars, however, the binary fraction increases to 
unity, as Sioo(Mi) m 1 and /m c = for Mi * c. The specific 
binary fraction then rises again to a maximum around Mi = d, 
and then rapidly drops to zero at Mi -2d- c. As a result of the 
rejection of low-mass companions, the overall binary fraction 
is always smaller than S. 

6.4.3. Pairing function SCP-III 

For SCP-III the mass ratio is redrawn when a companion with 
mass Mi < c is produced by the splitting algorithm (similar to 
PCP-III). This effectively corresponds to a (mass-specific) re- 
normalised mass ratio distribution ti q (q) in the range c/(Mc - 
c) < q < 1. The resulting overall mass ratio distribution for 
SCP-III is given by: 



/ ? ,all(<?) = k 



r-ld 

Jc(l + ] 



h'(q)f Mc (Mc)dMc , 



(24) 



(l+V?) 



where the lower integration limit is set by the condition M^ > c 
and k is again a normalisation constant. The expressions for the 
specific distribution f q -jt y (q) for a restricted primary mass range 
are identical to those for the SCP-II case (see Appendix I A. 4b . 
except that h' q (q) replaces h q {q) everywhere. As a result of the 
renormalisation of the mass ratio distribution, the mass ratio 
distribution is a function of spectral type. The specific mass 
ratio distribution for the lowest-mass cores is strongly peaked 
to q = 1 . 

The overall binary fraction S a n equals S for SCP-III. The 
specific binary fraction is given by Eq. ( 1221 . Note however, that 
the primary mass distribution /m,(M[) resulting from SCP-III 
is different from that of SCP-III; see Appendix IA.4.21 for de- 
tails. 

7. Differences between pairing functions 

As discussed in the previous sections, each pairing function re- 
sults in different properties of the binary population. In this sec- 
tion we provide an overview of the major differences and sim- 
ilarities. In general, each of the properties described below de- 
pends on the choice of the pairing function, a generating mass 
distribution /m(M), a generating mass ratio distribution f q (q), 



and a generating binary fraction S. Our example models have a 
Kroupa generating mass distribution with 0.02 < M < 20 M , 
a flat generating mass ratio distribution (if applicable), and a 
generating binary fraction of 100%. The resulting mass dis- 
tributions and mass ratio distributions for the different pairing 
functions are shown in Figs. [2] and [8] respectively. 

7.1. The mass distributions 

For pairing functions PCP-I, PCP-III and PCRP the primary 
mass distribution is identical to /m(M). For pairing function 
PCP-II, /m, (M) is very similar to /m(M), but it contains less 
primaries of low mass due to rejection of very low-mass com- 
panions (which mostly occurs among low-mass primaries). The 
latter "primaries" are considered as single stars after removal 
of their companions. For pairing function RP the primary mass 
distribution is more massive than /m(M) due to swapping of 
the components that are drawn from /m(M). Pairing functions 
SCP-I, SCP-II and SCP-III result in smaller average primary 
masses than /m(M) due to the the core splitting. 

For all pairing functions the companion mass distribution 
/m, (M) is shifted to lower masses with respect to the generating 
mass distribution / M (M). Pairing functions PCP-III, PCRP and 
SCP-III result in a large number of low mass companions; the 
companion mass distribution is sharply peaked at M2 ~ c; see 
Figs. [3] and [8] For a Kroupa generating mass distribution, the 
other pairing functions show a peak in the companion mass dis- 
tribution around M2 ~ 0.05 M . Pairing functions PCP-I and 
SCP-I may result in arbitrarily small companion masses, but all 
other pairing functions have Mi.min - c - All pairing functions 
have M2, max = d. For each pairing function described in this pa- 
per, the companion mass distribution depends strongly on the 
selected primary mass ra nge (see Fig. O, but is i ndepen dent 
of S. Note, however, that iMetchev & Hillenbrand! J2OO8) find 
evidence for a universal companion mass distribution among 
stellar and substellar primaries. 

We define the binary system mass distribution /m t (M) as 
the distribution of masses Mj — M\ + M2 of all binary systems 
in a population. For pairing function SCP-I the system mass 
distribution equals the generating mass distribution: /m t (M) = 
/m(M). For all other pairing functions the expression for the 
system mass distribution is different from /m(M), /m, (M) and 
/m, (M). For random pairing, for example, the system mass dis- 
tribution can be described as a convolution (Eq. [9]>. 

Pairing functions with & < 100% result in a population of 
single stars. For PCP-II and SCP-II, even those with £ = 100% 
result in single stars. In our analysis we discriminate between 
the single mass distribution /m s (M) and the primary/single 
mass distribution /m, s (M). The latter distribution includes both 
single stars and primaries and is important for observers. Target 
lists for binarity surveys are often defined by /m, s (M) for a cer- 
tain mass range, as in practice it is often unknown whether a 
surveyed star has a companion star. For most pairing functions 
the single mass distribution f Ms (M) is equal to the generating 
mass distribution f M (M), and the fraction of single stars among 
the systems is given by 1 - S, where S is the generating binary 
fraction. For PCP-II and SCP-II additional single stars are ere- 
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Fig. 8. The mass distributions resulting from the different pairing functions, for the models shown in Fig. [2] From top to bottom, 
the panels show the primary mass distribution /m, , companion mass distribution the system mass distribution /m t , the mass 
distribution of the primaries and single stars combined /m 1s (i.e., the targets in an observational binarity study), the single star 
mass distribution /m s , and the distribution of all stellar masses fM M ■ Each distribution is normalised so that its maximum is unity. 
The vertical axis has a linear scale (rather than a logarithmic scale; cf. Fig. [TJ so that the differences are clearly shown. Each 
column corresponds to a different pairing function, which is indicated at the top. 



ated due to the rejection of very low-mass companions. The 
generated single stars are mostly of low mass. Only for pairing 
functions PCP-I, PCP-II, PCP-III and PCRP, the primary/single 
star mass distribution /m, s (M) is equal to the generating mass 
distribution f M (M). The primary/single star mass distribution is 
biased to higher masses for RP, and to lower masses for SCP-I, 
SCP-II and SCP-III. 

The bottom row of Fig. [8] shows the mass distribution 
faii{M) of all stars, including primaries, companions, and sin- 
gle stars. This mass distribution, which includes all stars in the 
population, is referred to as the mass function of a stellar pop- 
ulation. For a zero-age population, f a ii(M) is an initial mass 
function (IMF), e.g., the Kroupa IMF (Eq. Only for pair- 
ing function RP is the individual star mass distribution f a ii(M) 



equal to the generating mass distribution /m(M); the other pair- 
ing functions result in a distribution f a ii(M) which is biased to 
lower values with respect to /m(M). 

7.2. The mass ratio distributions 

The resulting overall mass ratio distribution f q ^\\{q) depends on 
the pairing function. The top panels in Fig. [2] show the signif- 
icant difference between the overall mass ratio distribution for 
the different pairing functions. The overall mass ratio distribu- 
tion / 9j aii(<7) 1S equal to the generating mass ratio distribution 
fq{q) for pairing functions PCP-I and SCP-I. For pairing func- 
tions PCP-II, PCP-III, SCP-II and SCP-III the overall mass ra- 
tio distribution is biased to higher values of q with respect to 
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Table 2. The mass ratio distribution resulting from different 
pairing functions, for different samples, as compared to the 
generating mass ratio distribution h q (q). The columns show the 
pairing function, the overall mass ratio distribution, the specific 
mass ratio distribution for high-mass stars, the specific mass ra- 
tio distribution for low-mass stars. The symbols the table indi- 
cate whether the mass ratio distribution is equal to (=), almost 
equal to (»), or biased to low mass ratios ({) or high mass ra- 
tios (t ), with respect to h q (q). As the distribution h q (q) is unde- 
fined for RP and PCRP, the properties of / ?:J u [ (q) with respect 
to f q ,Mi(q) are indicated. 
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f q (q). For these pairing functions binaries with very low com- 
panion masses, and thus often very low mass ratios, are either 
rejected or redrawn, resulting in systematically higher values of 
q. The overall mass ratio distribution of pairing functions RP 
and PCRP are purely a result of the mass distribution /m(M) 
and depend strongly on its lower and upper limits. 

Only for pairing function PCP-I is the mass ratio distribu- 
tion independent of spectral type, and thus equal to the over- 
all and generating mass ratio distributions. For all other pairing 
functions the specific mass ratio distribution is a function of the 
primary mass. The middle and bottom panels in Fig.|2]show for 
each pairing function the specific mass ratio distribution for tar- 
get samples of different spectral types. For each of these pair- 
ing functions, high-mass binaries have on average a lower mass 
ratio than low-mass binaries. The lowest-mass binaries have a 
mass ratio distribution peaked to q = I. 

Figs. |2] and [6] illustrate the strong dependence of the spe- 
cific mass ratio distribution on the targeted sample in the sur- 
vey. Fig. |4] shows a generalised version of these figures. Each 
panel shows the two-dimensional distribution f(q,M\) for the 
different pairing functions. Care should thus be taken when ex- 
trapolating the results to the population as a whole. The in- 
terpretation of the observations is further complicated by the 
instrument bias and observational errors, an effect we will dis- 
cuss in §|9] An overview of the mass ratio distribution changes 
is presented in Table [2] 

7.3. The binary fractions 

The overall binary fraction T> d \\ is equal to the generating bi- 
nary fraction & for most pairing functions. Only for PCP-II and 
SCP-II the overall binary fraction is smaller than & because 
of the rejection of low-mass companions. Table [3] provides an 
overview of the changes in the binary fraction as a function of 
spectral type and primary mass range. Most pairing functions 



Table 3. The specific binary fraction as compared to 

the generating binary fraction S, for different the pairing func- 
tions. The columns show the pairing function, the overall bi- 
nary fraction S d \\, the specific binary fraction for high-mass 
stars SM,«M milx and for low-mass stars Sm^m,^,- The last col- 
umn shows whether the resulting binary fraction is equal to S 
for any primary mass range. See also Fig.|7]for several exam- 
ples. 
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result in a mass-dependent binary fraction. If the generating bi- 
nary fraction is smaller than 100%, all pairing functions except 
PCP-I, PCP-III and PCRP result in a specific binary fraction 
that depends on primary mass (see, e.g., Fig. [7). On the other 
hand, if £ = 100%, only pairing functions PCP-II and SCP-II 
result in a varying Sm, (Mi), due to the rejection of low-mass 
companions. 

For a sample of binaries with high-mass primaries, 
( Bm 1 (Mi) is approximately equal to £ for pairing functions 
PCP-I, PCP-II, PCP-III and PCRP. For RP, the specific bi- 
nary fraction for high-mass binaries is larger than S (unless 
S = 100%), while for SCP-I, SCP-II and SCP-III the binary 
fraction for high-mass stars is smaller than T> (unless £ = 100). 
For low-mass binaries, Sm 1 (Mi) is equal to T> for pairing func- 
tions PCP-I, PCP-III, and PCRP. For RP, the specific binary 
fraction for low-mass stars is only equal to the generating bi- 
nary fraction if the latter is & = 0% or S = 100%, and smaller 
in the other cases. For PCP-II the specific binary fraction for 
low-mass stars is smaller than S. For SCP-I and SCP-III the 
binary fraction for low-mass stars is larger than T>. For SCP-II 
the specific binary fraction for low-mass stars may be larger or 
smaller than S, depending on the properties of /m(M) and f q (q) 
and the value of S. 

8. Dependence on generating properties 

In this section we discuss how the properties of a binary pop- 
ulation depend on the attributes for the pairing functions: the 
generating mass distribution /m(M), the generating mass ratio 
distribution f q (q), and the generating binary fraction S. For the 
most important properties of a population (with respect to bina- 
rity), we list in Table |4] whether or not they depend on /m(M), 
f q (q), or < B, for each of the eight pairing functions described 
in this paper. Note that the system mass distribution /m t in- 
cludes both single stars and binary systems. The specific mass 
ratio distribution f q -M x (q) and specific binary fraction T>m, (Mi) 
in Table |4]only refer to samples where all binaries have identi- 
cal primary mass M\ . For a sample with a finite primary mass, 
range, the results for the specific mass ratio distribution and 
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Table 4. Do the properties of a population resulting from a certain pairing function depend on its generating properties? The three 
generating properties /m(M), f q (q) and S are listed in the top row. For each of the ten quantities listed in the left-most column, 
we list whether or not they depend on the choice of /m(M), f q (q) and S. We indicate dependence and independence with the 
symbols V an d x > respectively. The results for f q -M, (q) and Sm 1 (Mi) are valid only for a population of binaries with identical 
primary mass M\ . The distribution over total system mass fu T includes both binary systems and single stars. 



specific binary fraction are mostly identical to those of f q ^n(q) 
and Sail, respectively. Table [4] also illustrates which properties 
of the population can be used to recover /m(M), f q {q), and S. 
For example, for RP, the specific binary fraction T>m, (M\) pro- 
vides information on the generating mass distribution (e.g., the 
IMF). 

The generating mass distribution. For obvious reasons, all 
mass distributions listed in TableHJdepend on /m{M). The mass 
ratio distributions for RP and PCRP are defined by, and depend 
strongly on the properties of /m(M) and S. For PCP-I/II/III, 
fq;M, (q) does not depend on /m(M), as both Mi and q are drawn 
independently from their generating distributions. Obtaining 
the overall mass ratio distribution f qA \\{q) requires integration 
over the primary mass distribution; as for PCP-II and PCP-III 
fqMiiq) varies with Mi, so does f q ^\\(q). The overall binary 
fraction 2? a u is independent of /m(M), except for PCP-II and 
SCP-II, for which low-mass companions are rejected. The spe- 
cific binary fraction < Bm x (M\), does not depend on /m(M) for 
PCRP and PCP-I/II/III due to the independent drawing of Mi 
and M 2 (or q), while it does vary with /m(M) for SCP-I/II/III 
as Mi and q are not drawn independently, and for RP as Mi 
and M2 are not drawn independently (due to the swapping of 
the components; see § 16.11 ). 

The generating mass ratio distribution. The generating 
mass ratio distribution f q (q) is undefined for RP and PCRP. 
For the other pairing functions, the dependence of (q) and 
fg,aa{<l) ° n fg(q) is obvious. As companion masses are derived 
from q, the distributions /m 2 and /m m depend on the choice of 
f q {q). Note that that for SCP-II, all parameters vary with f q {q), 
and for PCP-II most parameters (except the primary/single 
mass distribution) vary with f q (q). The properties of the sin- 
gle stars do not depend on f q (q), except for PCP-II and SCP-II, 
where additional single stars are created due to the rejection of 
low-mass companions. 

The generating binary fraction. The dependence of S a ii and 

(Mi) on the generating binary fraction S is trivial. The dis- 
tributions that do not involve single stars, such as /^,, /m,, 



fq.au(l) an d fq;M, (<?)> by definition never depend on the choice 
of S. The mass distribution of all stars /m m (the "IMF") de- 
pends on the choice of £ for all pairing functions, except for 
RP. 

9. Interpretation of observations 

Binarity and multiplicity provide important information about 
the out come of the st ar forming process in different environ- 
ments ( Blaauw 196 lb . In this paper we explore this issue by 
making the assumption that binary stars are formed through a 
simple "pairing function". In reality the distribution of stars in 
/mi,Mo(Mi, Mz) is the result of complex physics involving the 
collapse of a molecular cloud into stars and stellar systems with 
disks (which can themselves fragment), followed by dynami- 
cal evolution of the protocluster (see §[5j . The resulting "pair- 
ing function" may thus not be describable in terms of the sim- 
ple probability distributions given in this paper (§|6j). However, 
even if we proceed from our assumption that the pairing of bi- 
nary stars involves the random selection of a primary mass fol- 
lowed by the secondary (RP, PCRP, PCP-I/II/III) or the random 
splitting of cores (SCP-I/II/III), the interpretation of the obser- 
vations is not trivial: 

- There is a large space of possible models for the forma- 
tion and evolution of a binary population. In this paper we 
describe eight pairing mechanisms. For each pairing func- 
tion there are a large number of possibilities for the gen- 
erating mass distributions /m(M), mass ratio distributions 
fq{q), and binary fractions S. In addition, we have no a pri- 
ori restriction on the plausible formation mechanisms (see, 
however, §[5)- 

- The observations of a binary population are generally lim- 
ited. The surveys are incomplete, and are affected by se- 
lection effects and observational biases, and often only a 
limited set of binary population parameters is measured. 

When the inverse problem of obtaining the binary formation 
mechanism from the data is so poorly constrained, it is not pos- 
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sible to find the "best model" for binary formation from a fit to 
the data. The only thing one can realistically do is to exclude 
models that are not capable of reproducing the data and accept 
that all other models offer plausible binary formation prescrip- 
tions. 



9.1. Constraints from observations 

Below we list several properties that have been identified for 
various binary populations over the last decades. These provide 
important information on the primordial pairing function, and 
the formation and evolution of binary populations. 

The observed mass-dependent binary fraction. The ob- 
served binary fraction is known to increase with i ncreas - 
ing primary s pectral t ype; se e , e.g. , Sterzik & Durisen j2004 : 
Kohler et alj d2006l) : lLadal d2006l) : iBouvetalJ feOOoF lor 
an overview. For early-type ( O/B/A) stars the binary frac- 
tion approaches 100% (e.g., Abt et alj Il990t iMason et al. " 



199 81 IShatsky & Tokovininl 120021: iKobulnickv & Frverj l2QQ7t 
Kouwenhoven et alj l2007bl) . The binar y fraction de c reases 

1971 



to 50 



60% for 



F/G-type 



stars (Abt & Lew 



Duq uennov & Mavorlll991|). For M-type stars the binary frac- 
tion is 30 - 40% (Fischer & Marcvlll992l: iLeinert et alll 19971 
Reid & Gizisl 11997^ and for late M-type stars and brown 



dwarf s the binary fraction decreases to 10 — 20% (Gizis et al. 



2003 



2003 



Close et al] 



Siegleretal 



2003 



2005 



Bouv et al. 120031: iBurgasser et al.. 
Maxted et alj|2008l) . Note that this 



correlation between mass and binary fraction is also pre- 
dicted by the hydrodynamical/sink particle simulations of Bate 
(2008). Assuming that this observational trend is not induced 
by selection effects, is inconsistent with pairing functions PCP- 
I, PCP-III and SCP-I, for which the binary fraction is indepen- 
dent of primary mass. Furthermore, observations have ruled 
out pairing functions RP and PCRP in various stellar group- 
ings (see §|5j. Pairing functions PCP-II, SCP-II and SCP-III 
remain options to describe the binary population, as for these 
pairing functions the binary fraction increases with increasing 
stellar mass, and the average mass ratio decreases with increas- 
ing mass. However, in this paper we merely describe simplis- 
tic (but frequently used) pairing functions. A deeper analysis, 
including a study of more complicated pairing functions, com- 
bined with further observations, is necessary for a full descrip- 
tion of the pairing function in the different stellar populations. 

Twin binaries. Observationally, there is a large prevalence 
of massive binaries with a mass ratio close to unity, often 
referred to as the "twin peak" in the mass ratio distribution 
dLucv & Riccolll979l: lTokovininll2000t IPinsonneault & Stanek 
20061: lLucyll2006l: ISoderhjelmll2007l) THigh-mass twin binaries 
are extremely rare for RP and PCRP. For pairing functions PCP 
and SCP, high-mass twin binaries only occur frequently when 
this is explicitly put into the generating mass ratio distribution. 
A high prevalence of low-mass twin binaries, on the other hand, 
naturally results from all pairing functions except for PCP-I. In 
general, peaks in the mass ratio distribution can occur for any 
mass ratio (see, e.g., Fig. [6]). The location of the peak depends 
on the pairing function and the primary mass range, and, if ap- 
plicable, the mass ratio distribution. In general, the peaks occur 



at low-<7 for a sample of high-mass binaries, and at high-g for a 
sample of low-mass stars. Pairing functions RP and PCRP are 
thus excluded, while pairing functions PCP and SCP can only 
result in massive twin binaries if the corresponding generating 
mass ratio distribution is strongly peaked to q — 1 . 

The brown dwarf desert. The brown dwarf desert is defined 
as a deficit (not necessarily a total absence) of brown dwarf 
companions, either relative to the frequency of stellar com- 
panions, or relative to the frequency of planetary companion s 
dMcCarthy & Zucker man 2004; G refher & Lineweaverl 2006). 
Theories have been developed that exp lain the existence of 
the b rown dwarf d esert using migration dArmitage & Bonnelll 
20021) or ejection dReipurth & Clarkd 1200 ll) of brown darfs. 
The latter scenario "embryo ejection" is most popular, and pre- 
dicts ejectino of brown dwarfs soon after their formation. In 
this scenario, brown dwarfs could be considered as failed stars. 
Kouwenhoven et alJd2007al) . however, show that the scarcity of 
brown dwarf companions among intermediate-mass stars can 
also be explained by an extrapolation of the mass ratio distri- 
bution into the brown dwarf regime; PCP-I/II/III are thus not 
excluded by the presence of the brown dwarf desert. 

The (initial) mass distribution. The initial or present-day 
mass distribution / a u (M) of a stellar population sets strong con- 
straints on the star formation process, and is an important fea- 
ture of each pairing function. The distribution i s often derived 
after its members are securely identified (e.g.. IKjoupe 12001 



Preibisch et ap2002l: lHarayama et ap200a IStolte et alJ2008h . 
The measured mass distribution is often the distribution of sin- 
gle/primary masses /m, s (M), as it is not known which mem- 
bers are single and which are binary, which results in a mea- 
sured mass distribution that is biased to higher masses with 
respect to the overall mass distribution f d \\{M), which, if mea- 
sured just after star formation, is the IMF. As stellar masses are 
often derived from measured luminosities, the presence of un- 
resolved binaries and crowding may further bias the measured 
(initial or present-day) mass distri bution to higher masses (e.g., 
Vanbeverenl 1981 iMafz Apellanizll2008l) . Over the last decade, 
considerable effort has been put into studyi ng possible envi- 
ronmental dep endences of the IMF (see, e.g., Elmegreenl2007 : 
Kroupal 20081 for an overview and examples). The IMF of a 
population is presumably a result of the form of the initial core 
mass function, and the primordial pairing fun ction (i.e., how 
these cores fragment into multiple systems, see lGoodwin et al ' 



120081) . An environmental dependence of the primordial pairing 
function (e.g., mass ratio distribution, binary fraction) implies 
a different outcome of the star formation process with environ- 
ment, it also almost certainly implies an environmental depen- 
dence of the IMF (unless the core mass function changes in 
such as way as to mask this change). 

9.2. Recovering the pairing function 

The pairing function, f M (M), f q (q) and & can in principle be 
derived from observations of binary systems, provided that the 
observations cover a large part of the parameter space {Mi, q); 
see §[8] A significant complication, however, is introduced by 
selection effects, in particular by detection limits that prevent 
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the detection of faint companion stars. As an example, a twin 
peak for high mass binaries would rule out all of the models 
presented in this paper if f q {q) is assumed to be fiat for PCP and 
SCP (see § 19. It . On the other hand, a twin peak at the low mass 
end only rules out PCP-I 0/the generating f g (q) is fiat). So how 
should we proceed when interpreting the observational data? If 
we assume for the moment that we somehow know that one of 
the eight pairing mechanisms discussed in this paper occurs in 
Nature, we can advise the following: 

- Make sure the survey of the binary population is complete, 
i.e., that all primary masses are sampled. As is clear from 
Fig. [2] looking only at the low mass stars does not allow 
differentiating PCP-II, PCP-III, RP, PCRP, and SCP-I/II/III 
(all show a twin peak). 

- Examine not only the overall distribution of a certain pa- 
rameter (such as q or the binary fraction) but also study 
how it varies with primary mass. Again, from Fig.|2]it can 
be seen that when considering only a single row of panels 
it is not possible to easily differentiate the pairing mecha- 
nisms. However, when looking at the overall and specific 
mass ratio distributions, the differences do become clear. 

- Examine the combined behaviour of each parameter (the 
mass ratio distribution and the binary fraction) as a function 
of primary mass. The combination of these parameters con- 
strains the possible pairing mechanisms significantly fur- 
ther. 

- List all mechanisms capable of reproducing the observa- 
tions as possible solutions to the inverse problem. Do not 
try to give a single answer if this is not warranted by the 
data. 

In reality the number of possible models is of course much 
larger, especially if we start from arbitrary probability distri- 
butions that are not constrained by an understanding of the 
physics of binary formation. This is illustrated by the simple 
example in § |9.3| which shows that allowing a power-law distri- 
bution for f q {q) means that only RP and PCRP can be excluded 
based on an observed fiat distribution of q. 

The only practical way of excluding models of binary for- 
mation based on observations is to treat the inverse problem 
with Monte Carlo methods where the observations are pre- 
dicted from the model and compared to the real observations. 
In this method sele ction and observational biases should be 
include d (see, e.g., Kouwenhoven! 20061 : Kouwenhoven et al 



2007bh . Starting from models based on probability distribu- 
tions for a set of parameters may not be the most useful way 
of constraining the formation mechanism for binaries as this 
leaves a lot of freedom. It is more fruitful (but also more diffi- 
cult) to start from actual physical models of binary formation 
and see if these are capable of reproducing the observations. 

A further complication occurs when one wants to recover 
the primordial pairing function, i.e., the pairing function that 
is present just after star formation, as the pairing function 
of a stellar population evolves over time as a result of sev- 
eral processes. During the first stages of star formation, the 
newly formed p roto-binaries ar e affected by pre-main sequence 
eigenevolution (Kroupa 1995a c.b]) due to interaction with the 
remaining gas in the circumbinary disk. Dynamical interactions 



Table 5. Suppose that a stellar population has an observed mass 
ratio distribution / 3 ,obs(<?) = 1 for qm m < q < 1. What is the 
pairing function and the generating mass distribution f q {q)l 
This table lists the exponent y q of the generating mass distribu- 
tion f q {q) °c q lc < that is most compatible with the observations. 
We assume that g m ; n = 0.1, and binaries with q < 0. 1 cannot be 
detected. We ignore the other selection effects. The numbers 
in the three columns represent the most compatible values of 
y q for the overall mass ratio distribution f q {q), and for the spe- 
cific mass ratio distribution of binaries with A/B primaries and 
brown dwarf primaries, respectively. 



Pairing function 


y q (all stars) 


y q (AB stars) 


7, (BDs) 


Observed 


0.00 


0.00 


0.00 


PCP-I 


0.00 


0.00 


0.00 


PCP-II 


0.35 


0.00 


1.30 


PCP-III 


0.55 


0.00 


1.80 


RP 


excluded 


excluded 


excluded 


PCRP 


excluded 


excluded 


excluded 


SCP-I 


0.00 


-0.40 


0.10 


SCP-II 


0.30 


-0.40 


1.50 


SCP-III 


0.55 


-0.40 


1.75 



can result in ionisation of binaries, the formation of new binary 
systems, and exchange interactions, and th us alters the pairing 
funct i on of a stellar population (see, e.g., Kroupa et al. 1999, 



Preibisch et al 



2003; Duchene et al. 200^ 



LID 

4; 



Kohler et al 



2001: 

120061: iReipurth et al.ll2007h . The pairing function also changes 
due to stellar evolution, which can change the mass of one or 
both of the compo nents of a binary system, and in some cases 
in a merger (e.g., ISills et allEooltlGaburov etai1 l2008). or in 



the ej ection of one o f the components during a supernova event 
(e.g., |Bla 



1961). The primordial pairing function can be 
constrained using t he technique of inverse dynamical popula- 
tion synthesis (e.g., Kroupa 1995albh . in which the outcome of 
Af-body simulations is compared with present-day binary pop- 
ulation in a Monte Carlo way. 

Finally we stress that any interpretation of observations of 
a binary population in terms of the formation of binaries should 
start by stating the assumptions one makes in order to restrict 
the number of solutions to explore. That is, in the context of 
what class of binary formation mechanisms are the observa- 
tions interpreted? 

9.3. An example - an observed flat mass ratio 
distribution 

Most pairing functions result in a mass ratio distribution that 
varies with primary spectral type. For this reason one has to be 
cautious when interpreting the observations of a sample of bi- 
naries. Given the observed dataset, what is the pairing function, 
and what is the generating mass ratio distribution f q {q)l The 
answer partially depends on the generating mass distribution, 
which we assume to be the Kroupa mass distribution for now. 
More importantly, the answer depends on the properties of the 
surveyed targets. In this example we analyse three cases: a sam- 
ple where all binaries are studied, a sample of spectral type A/B 
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targets, and a sample of brown dwarf targets. Using Monte- 
Carlo techniques we determine which pairing function is con- 
sistent with an observed flat mass ratio distribution / 9>0 bs(<7) = 1 
for each subset, and, if applicable, which generating mass ratio 
distribution. We assume that the generating mass ratio distribu- 
tion has the form f q (q) oc q y i . We further assume that the obser- 
vations are complete in the range q m i n < q < 1 with g m ; n = 0.1, 
and that no binaries with q < q^ are observed due to incom- 
pleteness. 

Table shows the best-fitting values of y q for each pairing 
function. For each for the three samples, pairing functions RP 
and PCRP are excluded with high confidence; these are unable 
to reproduce the observed flat mass ratio distribution. The best- 
fitting value for pairing function PCP-I is y q = for all sam- 
ples. This is not surprising, as for (q) - f q (q) for this pair- 
ing function. For PCP-II and PCP-III the derived y q for high- 
mass stars equals the observed value, but the other two samples 
contain more binaries with high mass ratios. For pairing func- 
tions SCP, the best-fitting intrinsic value of y q is smaller than 
the observed value of y q for high-mass binaries, but larger for 
low-mass binaries. This example illustrates that the inferred in- 
trinsic pairing properties may be significantly different from the 
observed pairing properties, depending on the pairing function 
and the selected sample of binaries. 

In practice, parameter distributions are often represented 
with a functional form. Suppose, for example, that we as- 
sume that the generating mass ratio distribution has the form 
f q (q) 00 q Jq - If we use this functional form for our model pop- 
ulation, and compare simulated observations with the true ob- 
servations for different values of y q , we will likely find a best- 
fitting y q . This does not necessarily mean that the generating 
mass ratio distribution has indeed the form f q (q) °c q 7tl . In this 
example we have added another assumption, i.e., that the mass 
ratio distribution has the form f q {q) °c q y '< . 

10. Summary and discussion 

We have described several methods of pairing individual stars 
into binary systems. We refer to these algorithms as pairing 
functions. Each pairing function is characterized by a gener- 
ating mass distribution /m(M) and a generating binary fraction 
S, and most additionally by a generating mass ratio distribution 
f q (q). Each pairing function results in a significantly different 
binary population. Depending on the pairing function and the 
mass range of the binaries studied, the resulting binary popu- 
lation may or may not have a mass ratio distribution or binary 
fraction that is equal to f q (q) or S, respectively. The binary 
fraction and mass ratio distribution generally depend strongly 
on the number of substellar objects in the population, and on 
the properties of the surveyed sample. 

Eight pairing mechanisms are discussed in detail. For ran- 
dom pairing (RP) both components are randomly drawn from 
the mass distribution /m{M). For primary-constrained random 
pairing (PCRP), both components are drawn from /m(M), with 
the constraint that the companion is less massive than the 
primary. For primary-constrained pairing (PCP-I, PCP-II, and 
PCP-III), the primary is drawn from /m(M), and the compan- 
ion mass is determined using a mass ratio distribution f q (q). For 



split-core pairing (SCP-I, SCP-II, and SCP-III), the core mass 
is drawn from /m (M), and the masses of the binary components 
are determined by the mass ratio f q {q), which splits up the core 
into two stars. The difference between the variants of pairing 
functions PCP and SCP lies in the treatment of low-mass com- 
panions (see §|6]l. Seven pairing functions naturally result in a 
specific mass ratio distribution that depends on primary spec- 
tral type, and five naturally result in a mass-dependent binary 
fraction. Seven out of eight pairing functions always produce 
a twin peak for low-mass binaries, while none result in a twin 
peak for high-mass binaries, unless the generating mass ratio 
distribution is strongly peaked to q — 1 . 

The differences between pairing functions are important for 
(i) the interpretation of observations, (ii) initial conditions of 
numerical simulations, and (iii) understanding the outcome of 
star formation: 

(i) The interpretation of observations The choice of the obser- 
vational sample may mislead the observer in deriving the over- 
all properties of a stellar population, as most pairing functions 
have a mass-dependent binary fraction and mass ratio distri- 
bution. If the binary fraction or mass ratio distribution of two 
samples (e.g., systems with B-type primaries and those with 
M-type primaries) are different, this does not necessarily mean 
that the underlying pairing function is different. A significant 
further complication is introduced by observational selection 
effects, which artificially decrease the binary fraction and in- 
crease the average mass ratio. The only practical way to ac- 
count for these is by using a Monte Carlo approach, and to 
compare simulated observations of a model population with the 
results of the binarity survey, taking into account all sampling 
and selection effects (§|9]l. 

(ii) Initial conditions for numerical simulations A choice for 
the pairing function has to be made when generating initial 
conditions for simulations of star cluster simulations with bi- 
naries. The simplest choice is random pairing (RP), although 
this pairing function is excluded from observations and not ex- 
pected from star formation. When modeling star clusters, one 
has to be aware that most pairing mechanisms result in mass- 
dependent properties, such as mass ratio distribution and binary 
fraction. The choice of the pairing function affects the outcome 
of the simulations, such as the dynamical evolution of star clus- 
ters, mass segregation, and the number of contact binaries and 
supernovae. 

(iii) Star formation. Different star formation scenarios result 
in different mechanisms of pairing stars into binary systems. 
After star formation, the pairing function is altered by dynam- 
ical interactions and stellar evolution. Random pairing, how- 
ever, is not predicted by star formation models, and is excluded 
by observations (see §[5j . The binary fraction and mass ratio 
distribution can be used to discriminate between the different 
pairing functions. Although the pairing functions described in 
this paper are common in literature, we do not suggest that one 
of these pairing functions indeed describes the natural outcome 
of the star forming process. It is, for example, possible that 
the pairing properties are a function of primary mass or core 
mass. If this is the case, it may indicate different formation pro- 
cesses for different masses. Due to the lack of observations, this 
has not been studied in detail, apart from the extreme ends of 
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the mass dist ribution (v ery massive stars and b rown dwarfs). 
Kroupa et al.l ( 12003b and lThies & Kroupal (12007b . for example, 
find that the observed IMF and binary population can be ex- 
plained by separate formation mechanisms for stars and brown 
dwarfs (see § 16. 1.1b . Nevertheless, the pairing functions de- 
scribed in this paper are useful tools to describe the outcome of 
star formation simulations. In order to constrain the primordial 
binary population from observations, one does not only have to 
take into account the selection effects, but also the change in the 
pairing function that has occured due to the effects of stellar and 
dynamical evolution. Over the last decade, considerable effort 
has been put into studying possible environmental dependences 
of the IMF. The IMF of a population is presumably a result of 
the form of the initial core mass function, and the primordial 
pairing functi on (i.e., how these c ores fragment into multiple 



systems, see iGoodwin et al.l 120081) . An environmental depen 



dence of the primordial pairing function (e.g., mass ratio distri- 
bution and binary fraction) implies a different outcome of the 
star formation process and almost certainly an environmental 
dependence of the IMF (unless the core mass function changes 
in such as way as to exactly mask this change); see §[8] 

Each pairing function, as well as each subset of stars, re- 
sults in a different mass ratio distribution and binary fraction. 
It is therefore of great importance to carefully study selection 
effects in observations, and to clearly state the pairing mecha- 
nism used in simulations, in order to make statements about the 
star formation process. The pairing functions described in this 
paper are likely too simplistic to describe a realistic stellar pop- 
ulation. However, they are frequently used to describe observa- 
tions and simulations. The next step forward is to fully charac- 
terize the binary population of several young stellar groupings; 
only in this way the star formation process can be recovered. 
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Appendix A: The pairing function and mass ratio 
distributions 

In this appendix we discuss for each of the pairing mechanisms 
described in the main paper how to calculate the mass ratio 
distribution function f q {q). The masses of primary (Mi) and 
secondary (M 2 ) are drawn from the same generating mass dis- 
tribution /m(M), or alternatively, the primary is drawn from 
/m(M) while the secondary is drawn from a generating mass 
ratio distribution. The generating mass distribution is treated 
throughout this appendix as a probability density: 



P(M <t) = F M (t) 



-f 



f M (x)dx , 



and 



f 



f M (M)dM = 1 



(A.l) 



(A.2) 



In the subsequent section we describe for each of the pair- 
ing functions first how to calculate f q (q) without specifying the 
generating mass distribution and then we work out the resulting 
mass-ratio distributions for the single power-law /m(M) which 
is given here as: 



f(M) = aM~ a with a 



yc' 



d l-a _ c \-a 



1 - (c/dy 



(A3) 



where a + 1 and y — a - I and c and d are the lower and upper 
limits on the mass distribution. 



A.1. Random pairing 

The simplest choice for a pairing function is that of 'random 
pairing' (RP). In this case both component masses are drawn 
independently from /m(M) and swapped, if necessary, to en- 
sure that M 2 < M\. In this case the joint distribution function 
frpiMi , Mi) for the component masses is given by: 



f lv (M l ,M 2 ) = 2/m,(M 1 )/ M2 (M 2 ) , 



(A.4) 



where the factor 2 accounts for the fact that the pairs of masses 
are swapped in order to ensure that Mi < Mi . The masses are 
restricted by a lower limit c and and upper limit d which leads 
to the domain of / lp (Mi, M 2 ) being defined as c < M\ < d and 
c <M 2 <M X . 

A.1 .1 . Mass ratio distributions for RP 

First, general expressions for the mass ratio distributions 
are derived before working out specific examples. To de- 
rive the mass ratio distribu tion we follow the appendix of 
Piskunov & MaFKovl (11991I) and derive f q {q) from its cumu- 
lative distribution function: 



P(y/x <q) = FJq) 



(x, y)dxdy 



ff f * 

s 

2/(x 

Jclq Jc 



)f(y)dydx , 



(A.5) 



where for ease of notation we use M\ = x and M 2 = y. The 
integration domain S is shown in the left-hand panel of Fig. [5] 
The probability density f q is then given by: 



(A.6) 



In observational surveys for binarity one is often restricted to 
a certain range of primary spectral types for observational rea- 
sons. To derive the effects of the selection on primary mass, the 
derivation of f q (q) again proceeds via the cumulative distribu- 
tion function, which is now given by: 



F q (q) = P(y/x < q\X] <x<x 2 ) 



Jfu, 



y)dxdy , (A.7) 



where the primary mass range is restricted to x\ < M\ < x 2 and 
the integration domain S' is as shown in Fig. [5] The integration 
domain limits now depend on whether q is larger or smaller 
than cjx\. For q < c/x\ the integration domain 5" is given by: 
c/q <x<x 2 Ac<y<qx (middle panel in Fig. |5), while for 
q > c/xi S' is defined by x\ <x<x 2 Ac<y<qx (right- 
hand panel in Fig . |5j ■ The minimum possible value of q is c/x 2 
in this case, implying that F(q) = for q < c/x 2 . Thus the 
expression for the mass ratio distribution now becomes: 



f 



F q (q) 



q < — 



kf c l£ qX 2f(x)f(y)dydx i < q < j- 
kf x X2 f c 9X 2f(x)f(y)dydx i<q<l 



(A.8) 



and the probability density f q is again derived according to 
equation dA.6l . The normalisation constant k can be calculated 
from the the condition F q (l) = 1. 

A.1 .2. Binary fractions for RP 

For RP, the overall binary fraction S d \\ is equal to the gener- 
ating binary fraction S. The specific binary fraction Sm,(Mi), 
however, is generally a function of primary mass. For a given 
primary mass, the number of single stars S (Mi) with mass Mi 
and the number of binary stars B(M\) with primary mass M\ is 
given by: 

S Mi = Sf M (M l )dM l , (A.9) 
where S is the total number of single stars in the system, and 



B Ml =2Bf M (M0dM 1 f M (M')dM', 



r 



(A. 10) 



where B is the total number of binary stars in the system and 
c is the minimum stellar mass. The specific binary fraction 
S Ml (Mi ) = B Ml /(S Ml + B Mi ) is then given by: 

^,(Mi)= 2BFM) J*L± + i)-\ (A.11) 
2BF M (Mi) + 5 \2Fu(Mi) ) 

where F M (M\) is the cumulative mass distribution (i.e., the 
primitive of f M ) evaluated at mass Mj. 

An example of the mass-dependent binary fraction result- 
ing from RP is shown in Fig. [7] Clearly, we have Smi(Mi) = 




Fig.A.l. Mass ratio distributions f q (q) for random pairing from the single power-law /m(M). The curves are shown for a = 
2.35, 1.30, 0.30, -2.50 and four ratios of the lower to the upper mass limit of f M (M): c/d = 0.01, 0.1, 0.3, 0.5. 



1 when 5=0, and S Ml (Mj) = when B = 0. Also, 
S Ml (M max ) = 2B/(2B + S) > £ and S Ml (M min ) = 0. The mass 
at which the specific binary fraction equals the generating bi- 
nary fraction can be found by solving S = SmSMi), which 
gives F m(M\) = 0.5. In other words, for pairing function RP 
the overall binary fraction can be found at the median stellar 
mass (M). For larger primary masses, & Ml (M\) is larger, and 
for smaller primary masses, & Mt (M\ ) is smaller. 

A.1 .3. Single power-law mass distribution 



From Eqs. ( IA.5b . ( IA.61 l. and ( IA.8I 1 mass ratio distributions for 
specific choices of the generating mass distribution can be de- 
rived. The single power-law defined in Eq.( IA.3b is considered 
here. Without restrictions on the primary mass range, the cu- 
mulative distribution for q is given by: 




■qx 

F q (q) = I I 2a 2 x- a y- a dydx 

'clq *J c 
,2 



a 



yl c 2y 



[(cldf y q- y + q y - 2(cld) y ) 



The probability density for q is then obtained as the derivative 
with respect to q of F q which after some algebraic manipulation 
leads to: 



f q (q) = 



y 



(i - iddyf 

For c/d <K 1 one obtains: 

(a - l)q°- 2 a > 1 



(q a - 2 - (c/d) 2y q- a ) . (A. 12) 



(1 - a)q a a < 1 



(A.13) 



These approximations are poor when a » 1 . 

Fig. lA.ll shows a number of examples of the resulting mass 
ratio distributions for different values of c/d and the power-law 
slope a. The values for a represent the Salpeter mas s distri- 
bution (2.35), the slopes at the lower mass end for the Kroupa 
(2001) mass distribution (1.3 and 0.3), and a possible slope at 
the very low mass end of the mass distribution, where the num- 
ber of stars increases with m. The latter value may occur in a 
multi-part power-law mass distribution with a real turnover at 
the low mass end. For a > 1 the mass distribution decreases 
with log m and for a < 1 it increases with log m. However this 
does not represent a real turnover in the mass distribution, the 




Fig. A.2. Mass ratio distributions f q (q) for random pairing from a single power-law /m(M). The solid curve shows the complete 
mass ratio distribution (for all binaries in the population). The other curves show what happens to the observed f q (q) if the 
primary mass M\ is restricted to x\ < M\ < xi- The curves are shown for a = 2.35, 1.30, 0.30, -2.50. The value of c/d is 0.01 
and x\ and X2 are listed in the panels. 



number of stars still increases as the mass goes down as long 
as a > 0. Note how the peak in the q distribution changes as 
the values of a and c/d are changed. For very small values of 
c/d one can see from the approximation (IA. 1 3b that f q will be 
fiat for a ~ and »»2 and that it will peak at low values of q 
for < a < 2. For values of a larger than 2 the number of low- 
mass stars is so dominant that high values of q are favoured 
(i.e., both Mi and M2 are likely to be small). Conversely, for 
a < the rise of the number of stars with m again favours high 
values of q. For < a < 2 the ratio of probabilities to obtain 
low or high-mass stars is such that drawing two equal mass star 
is unlikely thus favouring low values of q. 



while for c/x-i < q < c/x\ F q (q) is given by: 

Fqiq) - H j c qX x- a y- a dydx = ^ Hx- a (q l -x l - - c l -)dx 

K ± (±(«r^ 2 2r + c- 2 v) - c- V) ■ 

For cjx\ <q<\ the cumulative distribution for q is given by: 
Fq(q) K I I x~ a y~ a dydx 

*JX\ *J c 

X r 2 ^2^ y<KX ^ J ~ X ^ ~ cJ<<x i ~ x 7^j ■ 



To find the expression for f q when the primary mass range 
is restricted Eq. ( lA.8t has to be worked out for the single 
power-law mass distribution. For < q < c/x2 F q {q) - 0, 



The normalisation constant k for the probability density of q 
can now be found by substituting q = 1 in the last expression 
for F q (q): 

k - r 2 (^ 2 2y " *\ 2y ) " c-y(x? - xA ' (A. 14) 
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Now the expressions for f q (q) can be derived by taking the 
derivative of the integrals above: 



/,(?) = 



9<t 



(A.15) 



If x\ = X2 f q (q) will be proportional to q a for c/x\ < q < 1 
and zero otherwise: 



/,(?) = 







l-a 



l-(c/ Xl )> 



-q 



(A. 16) 



Fig. IA.2I shows examples of the behaviour of f q (q) when the 
primary mass range is restricted. The solid lines show the dis- 
tribution of q for a complete binary sample. The dotted lines 
show what happens if there is only a lower limit (> c) on the 
primary mass range. For a > 1 the high mass ratios are pref- 
erentially removed because low-mass primaries are removed. 
Conversely if there is only an upper limit on M\ (< d) only 
low mass ratios are removed and the resulting f q {q) is given 
by the dashed lines. The dot-dashed lines show a generic case 
with c < x\ < X2 < d. The latter case for very narrow primary 
mass ranges f q {q) will converge to f q {q) cc q~ a . For < a < 1 
the effect of mass selection is to remove the low values of q, 
thus flattening the distribution and moving the peak. For a < 
the mass selection does not have much effect. This figure illus- 
trates that the interpretation of mass ratio distributions in terms 
of random pairing is not straightforward unless the generating 
mass function is well known and the observations are indeed 
complete over a known primary mass range. 

The expression for the specific binary fraction Sm, (M\) re- 
sulting from random pairing is listed in Table IA.31 and, for 
the Salpeter generating mass distribution (a = 2.35), shown 
in Fig. [7] 

A.2. Primary-constrained random pairing 

In this case (PCRP) the primary and secondary are again drawn 
independently from the generating mass distribution, however 
for the secondary the condition M2 < M\ is imposed before 
drawing the secondary mass. That is the probability density 
fM 2 (M2) is re-normalised to the interval [c, M{\ (recall that c 
is the lower mass limit on the mass distribution). 

A.2.1 . Mass ratio distributions for PCRP 

Writing the re-normalised secondary mass distribution as 
/'m 2 (^2) = f'(y), the expression for the joint probability dis- 
tribution fpcrp(Mi, M 2 ) is: 

/pcrpCMj , M 2 ) = f Ml (Mi )/' M2 (M 2 ) , (A. 17) 

which is normalised to 1 due to the re-normalisation of 
fuSMi). 

To derive f q one can proceed as for the RP case. The inte- 
gration domain is as shown in the left-hand panel of Fig. [5]and 



the expression for F(q) can be written as: 

~d r*qx 



F(q) 



f f 

*Jcjq Jc 



f(x)f'{y)dydx. 



(A.18) 



Note that the normalisation constant for f'(y) depends on x. 
The expression for the restricted primary mass range becomes: 



f 



F q (q) 



q < — 



k£; q ^f(x)f'{y)dydx ±<q<j [ , (A.19) 
kf* f c qX f(x)f'(y)dydx ±<q<l 



where the integration domains for c/x2 < q < c/x\ and cjx\ < 
q < 1 are shown in the middle and right-hand panels of Fig. [5] 
respectively. 

A.2.2. Binary fractions for PCRP 

For PCRP, the overall binary fraction S a n equals S, and spe- 
cific binary fraction Sm,(Mi) is equal to the generating binary 
fraction S for any mass Mi . 

A. 2. 3. Single power-law mass distribution 

For the single power-law mass distribution, with lower and up- 
per mass limits c and d, the joint distribution of x = Mi and 
y - M2 is: 

/ P crp(*,y) = a\x~ a 'a 2 (x)y~ a , (A.20) 

where the normalisation constant for f(y) depends on the pri- 
mary mass x: 



a; = 



yc> 



1 - (c/d)y 



and a2(x) = 



yc' 



1 - (c/xy 



(A.21) 



A distinction has to be made between the cases a < 1 and a > 1 
as will become clear below. 
F q (q) is given by: 



<Jc(q *J c 

= fli I 

Jcl 



y 2 c 2y x -a 



(l - (cid)y)(\ - (c/x)y) 



■y a dydx 



q-?(c/xy - 1 



(c/ x y - 1 



dx . 



Using the substitution z — c/x (which implies dx = -(c/z 2 )dz) 
the integral can be written in a more convenient form and its 
solution can be written as linear combination of the terms z y 
and ln(l - z y ) for a > 1, while for a < 1 the terms z 7 and 
ln(z y - 1) are involved. 

The resulting expression for f q {q) is: 



f q (q) 



(A.22) 

This expression diverges as q — > 1 which is due to the low- 
mass end of the primary mass distribution. As Mi approaches 
c the values of q will increasingly all be close to 1. The rate 
at which f q diverges depends on the value of a. For a > 
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Mass ratios for PCRP and single power law IMF with slope: 2.35 



Mass ratios for PCRP and single power law IMF with slope: 1.30 




Mass ratios for PCRP and single power law IMF with slope: 0.30 



c/d = 0.010 
c/d = 0.100 
c/d = 0.300 
c/d = 0.500 




c/d = 0.010 
c/d = 0.100 
c/d = 0.300 
c/d = 0.500 




Mass ratios for PCRP and single power law IMF with slope: -2.50 




Fig.A.3. Mass ratio distributions f g (q) for PCRP from a single power-law mass distribution. The curves are shown for a = 
2.35, 1.30, 0.30, -2.50 and four ratios of the lower to the upper mass limit of the mass distribution: c/d = 0.01, 0.1, 0.3, 0.5. 



the mass distribution peaks at the low-mass end thus causing a 
rapid divergence. 

Fig. IA.3I shows four examples of the mass ratio distribu- 
tion for PCRP from a single power-law mass distribution. The 
curves are for the values of c/d of 0.01, 0.1, 0.3, and 0.5. 

For the restricted primary mass range Eq. ( IA. 19b has to be 
worked out. For c/x 2 < q < c/x\ the expression for F q is: 



r*i rq* y C r x - a 

jciQ Jc i - (c/xy 



->X2 f*qx 

c/q 

which can be worked out to: 



F q (q) ■ 



f 

Jc/q 



* 2 cVq-Vx- 2 ?- 1 



(c/x)y - 1 



-dx . 



This integral can be worked out in the same way as for the full 
primary mass range. For c/x\ < q < 1 the expression for F q 
becomes: 



P c 7 q- y x- 2 ?- 1 - x-y- 1 
F qW x , , — — ; dx , 

J jr. 



The expression for f q {q) for a > 1 is: 


/«(*)= {MH'W-MfJ-^ 

For a < 1 the expression for /„ is: 

' 

/«(*)= ^MW) + (i) 

The normalisation constant k is: 

k = 7 ((c/x l y-(c/x 2 yr 1 



q< 



I 7 ) %*9<5 ■ 

fen* *i 



(A.23) 



q < - 



-I 7 ) t<<l<T 



(A.24) 



(A.25) 



(c/xy - 1 



Again, if X\ = x 2 , f q (q) will be proportional to q~ a for 
c/x\ < q < 1 and zero otherwise and the expression is the 
same as Eq. dA. 161 >. 

Examples of what happens in the case of PCRP when a 
selection is done on primary mass are shown in Fig. IA.4I The 




Fig. A.4. Mass ratio distributions f q (q) for PCRP from a single power-law mass distribution. The solid curves show the complete 
mass ratio distribution (for all binaries in the population). The other curves show what happens to the observed f q (q) if the 
primary mass Mi is restricted to x\ < M\ < xi. The value of c/d is 0.01 and x\ and X2 are listed in the panels. 



behaviour is qualitatively the same as for the random pairing 
case. 

A. 3. Primary-constrained pairing 

For the pairing mechanism discussed now the assumption is 
that there is a physical process which sets the primary mass 
and the mass ratio of the binary, rather than setting the masses 
of primary and secondary. For the primary constrained pair- 
ing mechanism (PCP) the assumption is that M2 is determined 
from M\ through the mass ratio q. The probability densities for 
M\ and q are specified in this case and they are assumed to be 
independent. That is: 

f P c P (M u q) = f Mx {Mx)h q {q) , (A.26) 

where the generating mass ratio distribution is written as h q (q) 
in order to distinguish it from the observed mass ratio distribu- 
tion f q (q). The latter can be obtained by integrating f pcp (Mi,q) 
over M\ : 

/«(«) = f f P c P (M u q)dM l . (A.27) 



There are a number of choices one can make in generating 
a binary population from / pcp . The generating mass ratio distri- 
bution is assumed to be specified for the interval < q < 1 and 
c < M\ <d, which leads to the following three possibilities: 

PCP-I All values of q are allowed which means that for a given 
primary mass M\,Q < M2 < M\ . Thus binary systems with 
'sub-stellar' secondary components are also allowed. 

PCP-II All values of q are allowed but only binary systems for 
which the secondary is 'stellar' (i.e. M2 > c) are retained. 
This amounts to integrating / pcp over the range c/q < M\ < 
d and re-normalising the resulting distribution of q to 1 . 

PCP-III Only values of q for which M2 > c are allowed, i.e. 
c/M\ < q < 1. This is equivalent to re-normalising the 
generating distribution h q (q) to the interval [c/M\, 1]. 

A. 3.1 . Mass ratio distributions for PCP 

Again we first derive the general expressions for f q (q) before 
discussing specific examples. 
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Mass ratios for PCP-II, single power-law IMF: a- 2.35 




Mass ratios for PCP-II, single power-law IMF: a= 0.30 
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Fig. A.5. Mass ratio distributions f q {q) for PCP-II and III for a uniform generating mass ratio distribution and primary masses 
from the single power-law mass distribution for a = 2.35, 0.30, -2.50. The curves are shown for three ratios of the lower to the 
upper mass limit of the mass distribution: c/d = 0.1, 0.3, 0.5. 
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PCP-I In this case one always obtains f q (q) = h q (q) because 
the distributions of q and M\ are independent. Also when re- 
stricting the primary mass range the observed mass ratio distri- 
bution is equal to the generating distribution. 

PCP-II Here the systems with M2 < c are discarded and then: 



h = k f / 

Jc/q 



h q (q)f Ml (M{)dM x , 



fq(q) = k f fp Cp (Mi,q)dMi 

Jc/q 

(A.28) 

where A: is a normalisation constant which ensures that 
J f q (q)dq — 1. Of course f q (q) = for < q < c/d. For a 
restricted primary mass range, c < x\ < M\ < X2 < d, the 
expression for f q is: 



f q (q) = 



0<q< j; 

k Qh q (q)f Ml (M 1 )dM l i_<q<j^ 
k h q (q)f Ml (M 1 )dM 1 f<q<l 



(A.29) 



where k is again a normalisation constant. For q < c/.X2 the 
value of c/q is larger than X2 so the lower integration limit of 
the integral in Eq. (1A.28E becomes X2, i.e., all companions with 
q < c/x2 are rejected, as their mass is smaller than the min- 
imum mass c. For q > c/xi we have c/q < x\ so the lower 
integration limit of the integral should be fixed at x\ . 

PCP-II I In this case the generating mass ratio distribution is 
re-normalised to the interval \cjM\ , 1 ] , resulting in a generating 
distribution h! q (q). The expressions for f q are then derived as 
for the PCP-II case: 



Jc/c 



f pcp (M u q)dM i 



Jc/q 



,(q)f Ml (Mi)dMi , 



(A.30) 

where a normalisation constant is now not needed (h' q (q) is 
normalised). For the restricted primary mass range the expres- 
sion is: 



f q (q) 



0<q<j; 
k^ q h' q {q)f Ml (M l )dM l j;<q<f x 
k^h' q (q)f Ml (M{)dMi X < q < 1 



(A.31) 



A.3.2. Binary fractions for PCP 

PCP-I The overall and specific binary fractions are always 
equal to the generating binary fraction S for PCP-I 



PCP-II Due to the rejection of low-mass companions with 
M2 < c, i.e., with q < q m i n (Mi) = cjM\, the binary fraction 
varies with primary mass: 



Sm,(Mi) 



h q (q) dq < S , 



(A.32) 



4mi„(A?l) 



where h q (q) is the generating mass ratio distribution. Note that 
HmSMx) is independent of the generating mass distribution. 
An example of the mass-dependent binary fraction resulting 



from PCP-II is shown in Figure [7] For high-mass binaries with 
Mi « d, very few companions are rejected as ^ m i n (Mi) <K 1, 
and hence Sm, {Mi) sb < B. For the lowest mass binaries in 
the sample, on the other hand, g m j n (Mi) m 1, and therefore 
Sm, (Mi) ~ 0. The overall binary fraction can be found by in- 
tegrating over primary mass 

S a u = J &M i (Mx)fM i (Mx)dM x , (A.33) 
and is always smaller than S. 

PCP-I 1 1 The overall and specific binary fractions are always 
equal to the generating binary fraction S for PCP-III 

A. 3. 3. Uniform mass ratio distribution and single 
power-law mass distribution 

The following is assumed for h q and /m,(Mi) = f(x) (using 
x — Mi for ease of notation): 

h q (q) = 1 < q < 1 
f(x) — ax~ a c < x < d , 

where a = yc y /(l - (c/d) 7 ) and a + 1. 

The PCP-II case For the full primary mass range we have 
from Eq. (IA281 : 



C d _ q 1 
fq(q) — k I ax a dx — k— 

Jc/q 1 



7 - (c/dy 



(c/d)y 



(A.34) 



c 1 

The expression for k can be found by solving J, d f q (q)dq = 1 
for k which leads to the final expression for the observed mass 
ratio distribution: 



fq(q) = 



< q < § 

(y+\)(qy-(c/dY) £ < „ < 1 

l-(y+l)(c/dy+y(c/dy+> d ~ " ~ 



(A.35) 



If the primary mass range is restricted to X\ < q < X2 the 
expression for f q can be derived using Eq. ( |A.29b : 

0<q<j; 
f q (q) = • Kq y ~ (c/x 2 y) i^q<i , (A.36) 
,*((c/xiF - (c/x 2 y) ±<q<l 

where the normalisation constant is: 

k= 111 

(y + 1) [(c/xi)r - (c/x 2 )r] - y [(c/x^ - (c/x 2 )r +l ] ' 

(A.37) 

Note that if xi — X2 the distribution of q will be uniform on 
the interval [c/x\, 1], i.e. f q (q) = 1/(1 - c/xi). The equations 
above do not hold for a = 0. This case is easily evaluated and 
the result for a = is: 



fq(q) = 



1 



1 - c/d + (c/d) \n(c/d) 



1 - 



c/d 



(A.38) 
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Fig. A.6. Mass ratio distributions f q (q) for PCP-II and III for a uniform generating mass ratio distribution and primary masses 
from a single power-law mass distribution. The solid curve shows the complete mass ratio distribution (for all binaries in the 
population) for c = 0.2 and d = 2.0. The other curves show what happens to the observed f q (q) if the primary mass M\ is 
restricted to x\ < M\ < x 2 . 
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The PCP-III case Now the mass ratio distribution is re- 
normalised to the interval [c/x, 1] which gives: 



h' q {q) = - 



1 



c/x 



(A.39) 



For the full primary mass range f q is obtained from Eq. (lA.30b 
[0 0<q<t 

fq(<l) = \ (J _t- ■ 

\ a l lq — dx 5^<7<1 
Restricting the primary mass range leads to: 
< q < ± 

k f* 2 —dx < q < 

Jcjq X-c -V.. ^ X\ , 



fM 



(A.41) 



ox 



-dx < q < 1 

Ci X-C Xt i 



The expressions for the PCP-III case can be written in a slightly 
more convenient form when using the substitution z — c/x. For 
the full primary mass range: 



/,(«) = 



< q < § 



(A.42) 



and for the restricted range: 



/,(?) 



1 X2 

, H — < q < — 

l-(c/d)? Jc/x 2 l-Z x 2 * X] 
e/.ti 

l-fc/d)? Jc/A'2 



(A.43) 



The integrals involving the term z y ~' /(l - z) can be solved and 
the general expression involves the hypergeometric function 
2-Fi. The expression is: 



1-z y y+l 



From the properties of th e hypergeometric function (see 
Gradshtevn & Rvzhikll2007l Chapt. 9) it follows that because 
(y+l) + l — (y+2) = Othe expression above converges through- 
out the unit circle in the complex plane except at \z\ = 1 . So 
there will be a singularity at q = 1 for the full primary mass 
range case. 

For integer values of y (or a) special care should be taken. 
From the expression above it is clear that y — or y — -1 
should be treated separately. Furthermore the hypergeometric 
series for 2F\(a,b;c;z) is indeterminate for c = —n where 
n = 0, 1,2, • • ■ if neither a nor b is equal to -m (where m < n 
and m is a natural number). Here this means that all cases 
y + 2 — -n should be treated separately which combined with 
the condition y + 0,-1 implies that whenever y = — n (i.e. 
a = -n + 1) the expression above will not apply. The case 
a — 1, y = 0, is excluded, so the special cases are y = -n - 1, 
a = -71, and then the expression for the integrals becomes: 



dz 



z n+ \\ -z) 



= y ± 

<H (n + 1 - Jfc 



i<=\ 



k)z : 



n+l-k 



-In 



1 -z 



(A.45) 



where the solution can be found from formula 2.117(4) in 
Gradshtevn & Rvzh ik (2007). In the integrals above the inte- 
gration constants were left out. 

For restricted primary mass ranges the results above can be 
used to evaluate the integrals listed in Eq. ( IA.43I ). The normal- 
isation constant k can be obtained from numerical integration 
of f q {q). Again, if X\ = X2 one obtains f q {q) = 1/(1 - c/x{) 
for q in the interval [c/x\, 1]. Figs. |A.5l and |A.6| show examples 



(A.40) of the resulting mass ratio distributions for h q {q) = 1 and vari- 
ous single power-law mass distributions for the primaries. The 
hypergeometric fu nction was calculated using the routine from 
Press etal.ld2007i §6.13). 

Finally, we have listed the expressions for the binary frac- 
tion as a function of primary mass in Table lA.3l 



A.4. Split-core pairing 

The mechanism of split-core pairing (SCP) works on the as- 
sumption that binaries are formed by the splitting of star- 
forming cores into two components. The component masses 
are specified through their mass ratio. The distribution func- 
tion for the core masses M c is given by the core mass function 
fc(M c ) and the mass ratio distribution f q (q) is specified inde- 
pendently. The masses of the primary and secondary are then 
given by: 

Mi = —^—M c and M 2 = — — -M c , (A.46) 



1 + q ~ 1 + q- 1 

where for simplicity the star formation efficiency is assumed to 
be 100% once the core mass is set. For a constant efficiency as 
a function of core mass this assumption has no influence on the 
results. Furthermore it is assumed that the minimum core mass 
is large enough to ensure that the primary is always of 'stellar' 
mass. That is the minimum core mass has to be at least twice 
the minimum stellar mass. 

The joint probability density / scp for M c and q is written as: 

/sc P (M c , q) = MM c )h q (q) , (A.47) 

where again the generating mass ratio distribution is written 
as h q (q) in order to distinguish it from the observed mass ra- 
tio distribution f q (q). The latter can be obtained by integrating 

/ S c P (M c ,^)overM c : 

/«(«) = J fsop(M c ,q)dM c . (A.48) 

As for PCP there are three cases: The generating mass ratio 
distribution is assumed to be specified for the interval < q < 1 
and we assume 2c < M c < 2d (with c and 2d being the mini- 
mum and maximum stellar mass, respectively), which leads to 
the following three possibilities: 

SCP-I All values of q are allowed, so that for a given core 
mass M c , Mi = M c /(1 + q) and < M 2 < Mi. Thus binary 
systems with 'sub-stellar' companions are allowed. 

SCP-II All values of q are allowed, but only binaries for which 
the secondary is 'stellar' (i.e. M 2 > c) are retained. This 
amounts to integrating / scp over the range c( 1 + 1 /q) < M c < 
2d and then re-normalising the resulting distribution of q to 
1. 
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SCP-III Only values of q for which M2 > c are allowed, i.e. 
cj (M c -c) < q < 1 . This is equivalent to re-normalising the 
generating distribution h q (q) to the interval [c/(M c - c), 1]. 

If we would have chosen our core mass distribution in the in- 
terval c < M c < 2d rather than 2c < M c < 2d, we would have 
encountered a further complication. For pairing function SCP- 
II the resulting single star mass is then occasionally smaller 
than c, and for SCP-III the splitting up is not possible if the 
binary system mass is smaller than 2c (see §|6]for details). For 
simplicity we avoid this issue in our analysis below, and sim- 
ply draw masses from the core mass distribution in the range 
2c<M c < 2d. 

A. 4.1 . Mass ratio distributions for SCP 

We give the general expressions for f q (q) before discussing 
specific examples. For ease of notation we use x = M c and 
fix) = / C (M C ). 

SCP-I If no selection on primary mass is made one always 
obtains f q (q) = h q {q) because the distributions of q and M c are 
independent. Unlike PCP-I, restricting the primary mass range 
now does cause changes of the observed mass ratio distribu- 
tion with respect to the generating one. Several cases have to 
be distinguished on the basis of the value of the primary mass 
selection limits x\ and X2 compared to the values of c, d, 2c, 
and 2d. Note that x\ and xi are limits on M\, not M c . 

In all cases discussed below a selection x\ < M\ < X2 in 
principle translates to X\(l + q) < x < X2(l + q) but the upper 
and/or lower limits on x and q used for integrating / scp and 
normalising f q (q) are different for each case. We introduce the 
following variables to distinguish the different cases: 

qo = 2c lx 2 - 1 q2 - 2d 1x2 - 1 

q\ — 2cjx\ - 1 qi — 2d/x\ — 1 

The value of x (core-mass) is restricted to [2c, 2d], hence for 
q < qo or q > qi the probability distribution for q vanishes (as 
these conditions imply M\ > X2 or Mi < x\, both of which are 
not allowed). For q < q\ the lower integration limit for x when 
determining f q is fixed at 2c and for q > q2 the upper limit is 
fixed at 2d. 

Case \: x\ < c < X2 < 2c In this case the value of qo is 
between and 1 and all other values of q, are larger than 1 . This 
means that two cases should be distinguished for f q , q < qo and 
q > qo, and that the lower integration limit for x is always fixed 
at 2c. The expression for f q is then given by: 



expression for /„. The lower integration limit for x is still 2c: 







< q < q 



fq(@) ~ I rx 2 (l+q) 

I k he h q (q)f(x)dx qo<q<\ 

where A: is a normalisation constant which follows from: 
-1 



(A.49) 



J 



f q (q)dq = 1 



(A.50) 



Case 2: xi < c A 2c < X2 < d This is the same as case 
1 except that now go < so that there is only one part to the 



fa® = k 



Jlc 



h q (q)f(x)dx < q < 1 



(A.51) 



where the normalisation constant k is obtained by integrating 
/.over [0,1]. 

Case 3: xi < c A d < X2 < 2d Now the value of q2 
is in the interval [0, 1] which means that for the case q > q2 
the upper integration limit for x is fixed at 2d which leads to a 
plateau in the probability density f q : 



f q (q) 



k Sic q h q(q)f(x)dx < q < q 2 

rid 

k hc h q (q)f(x)dx q 2 < q < 1 



(A.52) 



where the normalisation constant k is obtained through integra- 
tion over the intervals [0, ^2] and [^2,1]- 

Case 4: c < x\ < X2 < 2c Now the values of qo and q\ 
are in the interval [0, 1] which means that the lower integration 
limit for x depends on q. The probability density for q now 
consists of three parts: 



/,(?) = 







< q < qo 



kj 2 ~ h q (q)f(x)dx qo < q < q\ 
k f**^ hq(q)f(x)dx q\<q<\ 



(A.53) 



where k now has to be obtained from condition ( lA.50b by inte- 
grating over the intervals [qo,q\\ and [q\, 1]. 

Case 5: c < X] < 2c A 2c < X2 < d The value of q Q 
becomes less than zero and we have: 



fM) = 



kflc q h q(<i)f( x )dx < q < qi 



X \llq) ll l( C l)f( X ) dx <?1 < <? < 1 



(A.54) 



where k is now has to be obtained from condition ( IA.50t by 
integrating over the intervals [0, ^i] and [q\, 1]. 

Case 6: c < xi < 2c A d < X2 < 2d Now the values of q\ 
and q2 are in the interval [0, 1]. For q < q\ the lower integration 
limit is fixed to 2c and for q > q2 the upper integration limit is 
fixed to 2d. In addition it can happen that q\ < q2 or q\ > 
q2- The latter case will cause a plateau of constant probability 
density f q for q\ < q < q2- So now there are two 'sub-cases'. 

Case 6a: q\ < q2 



f q (q) 



k£f +g h q (q)f(x)dx < q < q\ 

k f*^+ q q J h q (q)f(x)dx qi < q < q 2 , (A.55) 

rid 

k 4(t+ ? ) hq(q)f(x)dx q 2 <q<l 



where the normalisation constant k is obtained through integra- 
tion over three intervals [0, q\], [#1, #2], an d [q2, 1]- 
Case 6b: q\ > q2 



f q (q) 



-x 2 (l+q) 
2c 
rid 

k Jlc h q(q)f( x )dx qi<q <q 

Id 

[(!+?) 



k *' h q (q)f(x)dx < q < q 2 

Id 
1c 

rid 

k Jx,n+^ h q (q)f(x)dx q 2 < q < 1 



(A.56) 
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where the normalisation constant k is obtained through integra- 
tion over three intervals [0, q 2 ], qi], and [q it 1]. 

Case 7: 2c < x\ < x 2 < d Now the values of go and q\ 
are less than zero and both q 2 and qi are larger than 1 . Hence 
the probability density consists of one part only: 



f q (q) = k 



I h q (q)f(x)dx < q < 1 . 

JxdU 



(A.57) 



1+3) 



where k is obtained by integrating f q over [0,1]. 

Case 8: 2c < x\ < d A d < X2 < 2d The value of q2 is 
now in [0, 1] so f q will consist of two parts: 



I k $Zl+l h q(q)f( x ) dx 0<q<q 2 

{ k Jud+g) h q (q)f(x)dx qi<q<i 



(A.58) 



where k is obtained by integrating f q over [0, q 2 ] and [q 2 , 1]. 

Case 9: d < x\ < x 2 < d Now the values of both q 2 and 
q-i are in [0, 1] and the probability density will vanish if q > qy. 



/«(«) 



k £i~(l+q) h q(q)f( x ) dx 0<q<q 2 
k J^ 1+q) h q {q)f(x)dx q 2 < q < <?3 ' (A.59) 
qi<q<\ 
where k is obtained by integrating f q over [0, q 2 ] and [q 2 , q^]. 

SCP-II The expression for f q is obtained by integrating the 
joint distribution / scp over x over the range c(l + 1/q) < x < 2d 
and normalising the resulting expression for f q to 1 : 



f q (q) 



JrOd 
e(l + ] 



(1 + 1/9) 



h q (q)f(x)dx , 



where k is obtained from the condition: 



f 

Jc (2 



fq(q)dq = 1 , 



(A.60) 



(A.61) 



c/(2d-c) 



with c/(2d - c) being the minimum possible value for q. 

When the mass range is restricted there are again a num- 
ber of cases to consider, depending on the values of x\ and 
x 2 . However the situation is less complicated than for SCP-I. 
First of all the value of c(l + 1/q) is always larger than 2c 
for < q < 1 which means that the value of q\ plays no 
role. Secondly the minimum possible value of q is c/x 2 and 
this quantity is always larger than qo for x 2 > c (which is 
mandatory) and therefore also the value of qo plays no role. 
For q < c/x 2 f q (q) is always zero. The values of q 2 and ^3 do 
matter as discussed below. 

Case 1: xi < x 2 < d In this case the integration limits for 
x, x\(l + q) and x 2 (l + q) are guaranteed to be less than 2d. For 
q < cjx\ the value of xi(l + q) is less than c(l + l/q) so the 
lower limit for integral over x is then fixed at c(l + l/q). The 
expression for f q {q) becomes: 



f 



/«(«) = 



< q < 



k laSl! h M)f( x )d x i<q<i , (A.62) 



k C^ h ^)mdx s.<q<\ 



where k is obtained from the condition: 

fq(q)dq = 1 ■ 



f 

*J cjxi 



(A.63) 



Case 2: x\ < d A x 2 > d Now the value of q 2 is less 
than 1 and for q > q 2 the upper limit of the integral over x is 
fixed at 2d (the value of x 2 (l + q) being larger than 2d). Now, 
q 2 > c/x 2 but q 2 may be larger or smaller than c/x\. So there 
are two sub-cases: 

Case 2a q 2 < c/x\ 



fq(q) 





Jc(l + 



0<q<j- 



*2<l+«) , , s r, S , C . 

H + l/ q ) h 9Wf( X ) dx T 2 ^1 < 12 

rid 

k Mnif q ) h 9(9)f(x)dx q 2 <q<± 

rid 



{AM) 



where k is obtained by integrating f q over [c/x 2 , q 2 ], [q 2 , c/xj, 
and [c/xi, 1], and applying condition (IA.63t . 
Case 2bq 2 > 



fq(q) 







0<^< ^ 

1 Xn 



k G) h ^)mdx f 2 <q<f i 

k Ix'o+q 9 ) h q(q)f( x *) dx T x ^q<qi 

rid 

k Jx,(i+q) h q(qW x ) dx qi<q<i 



(A.65) 



where k is obtained by integrating f q over [c/x 2 ,c/x\], 
[c/x\,q 2 ], and [q 2 , 1], and applying condition ( IA.63b . 

Case 3: x t > d This is the same as case 2 except that 
#3 < 1 which means that f q (q) = for q > q^ (this is due to 
the value of xi(l + q) becoming larger than 2d). There are the 
same two sub-cases: 

Case 3a: q 2 < c/x\ 



fq(q) 



*2(l+9) 
(1 + 1/9) "« 



0<q<j- 
h q (q)f(x)dx j-<q <q 2 



f 

*JE 
k f 2d 

rid 

k Uu q) h MW dx 







h q (q)f(x)dx q 2 <q < 

i<q<qi 

q-j < q < 1 



(A.66) 



where k is obtained by integrating f q over [c/x 2 , q 2 ], [q 2 , c/xj, 
and [c/xi, q$], and applying condition ( IA.63b . 
Case 3b q 2 > c/x\ 



0<q< 



fM 



k £l+9) h Mf(x)dx jr [ <q<q 2 
rid 

k J x h q (q)f(x)dx q 2 <q<q-i 



(A.67) 







qi < q < 1 



where k is obtained by integrating f q over [c/x2,c/xi], 
[c/xu q 2 ], and [q 2 , q^], and applying condition (IA.63t . 
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S' 




SCP-I The number of binary systems with a mass M\ is given 
by NB/m, (Mi )dM\ , where N is the total number of systems 
(singles plus binaries) in the population. The number of single 
stars with a mass M\ is given by N(l - B)f Mc (Mi)dMi. The 
specific binary fraction is thus: 



M 1 

Fig.A.8. The integration domains for obtaining the primary 
mass distribution and specific binary fraction for pairing func- 
tions SCP-I, SCP-II and SCP-III. The domain S for SCP-I is 
the parallelogram enclosed by the solid lines. The domain 5" 
for SCP-II and SCP-III is indicated with the shaded region. 



SCP-III Now the generating mass ratio distribution h q is re- 
stricted to the range c/(M c - c) < q < 1 and re-normalised. 
The corresponding distribution is h! q (q) and the expression for 
f q becomes: 



fM 



h> 

JcQ.+1/q) 



(q)f(x)dx, 



(A.68) 



where the lower integration limit is set by the condition M^ > c 
and the distribution is normalised. 



When the primary mass range is restricted to x\ 



< x < X2 



the cases and expressions for f q (q) are the same as for the SCP- 
II case, except that h'Jq) replaces h q (q) everywhere. 

A.4.2. Binary fractions for SCP 

In order to find the specific binary fraction as a function of pri- 
mary mass, we first need to find the primary mass distribution 
/m, (M\), which can be calculated as: 



/-?max(Ml) 

/m,(Mi)= f Muq (Muq)dq, 



(A.69) 



where q m i n (Mi) and <7 ma x(^i) are the minimum and maximum 
mass ratios for a given primary mass M\, and fy^AMi, q) is 
the joint probability density function for M\ and q. The lat- 
ter can be derived from the generating (core) mass distribution 
fM c (Mc) and the generating mass distribution h q (q): 



f Muq (M u q) = f Mc (M c (Mi,q))h q (q){\ + q) , 



(A.70) 



S Ml + B Ml B/m, (M\) + (1 - B)f Mc (M l ) 

(A.71) 

The integration limits for /m, (Mi) are indicated with domain S 
in Fig. IA.8I and are given by: 



[2- Mi/c c < Mi < 2c 







2c < Mi < 2d 



and 



1 c <Mi <d 

h- Mxld d<Mi<2d 



(A.72) 



(A.73) 



An example for the specific binary fraction resulting from SCP- 
I is shown in Fig. [7] The specific binary fraction for c < Mi < 
2c equals unity, as fu c (M{) — in this mass range. In the mass 
range 2c < Mi < d the specific binary fraction is practically in- 
dependent of Mi . Beyond Mi = d, the specific binary fraction 
rapidly drops to zero at Mi - 2d. The overall binary fraction 
Sail is always equal to B for SCP-I. 



SCP-II The integration domain for SCP-II to obtain fM^Mi) 
is indicated with region S ' in Fig. IA.8l i.e., for c < Mi <2d-c. 
The integration limits q m i n and q max are given by: 



and 



q mia = c/M\ c < Mi < 2d - c , 



[ 1 c < Mi < d 

\2- Mi/d d < Mi <2d-c 



(A.74) 



(A.75) 



The specific binary fraction for a population with a generating 
binary fraction of unity, i.e., B = 100%), is given by: 



»ioo(M,) 



r 



Cm f M ^ {Muq)dq 



(A.76) 



where the joint probability distribution is given by Eq. ( |A.70| i. 
For a given primary mass M\, the value of Sioo(Mi) given by 
the ratio between the (weighted) lengths of the horizontal line 
segments of domain S ' and S in Fig. IA.8l The integration limits 
<?min,2 and q m ax,2 are thus given by Eqs. (lA.74b and ( IA.75b . and 
the limits g m i n ,i and g ma x,i are given by Eqs. (IA.721 > and ( IA.73b . 
respectively. For systems with an arbitrary value of B, the num- 
ber of binary systems is given by 



B Ml =NBB m (Mi)f Ml (Mi)dMi, 
and the number of single stars by 



(A.77) 



where the factor (1 + q) is the Jacobian of the transformation S M{ = N(l -B)f Mc (Mi)dMi +NB(l -Bioo(Mi))f Ml (Mi)dMi , 
Mi = M c (q + I)" 1 and q = q. (A.78) 
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where the first term refers to the fraction of clumps that were 
assigned to be single, and the second term refers to the pri- 
maries that have become single due to the rejection of their low- 
mass companions. The specific binary fraction is then given by: 



B 



M | 



SiooCMOS/m^MO 



1 Mi + B Ml Sf Ml (Mi) + (1 - S)/m c (M,) ' 

(A.79) 

An example of the specific binary fraction resulting from SCP- 
II is shown in Fig. [7] For c < M\ < 2c, the above equation 
reduces to S(Mi) = $ m (Mi) as f Uc {M\) = 0. The binary 
fraction among the lowest mass stars (Mi = c) equals unity, 
after which it drops until M\ = 2c. Beyond that minimum mass 
it rises again and reaches its maximum at M\ = d, after which 
it drops again go zero at M\ - 2d, for the same reasons as for 
the SCP-I case. The overall binary fraction can be calculated by 
integrating T>(M\ ) over primary mass, weighed by the primary 
mass distribution /m,(Mi), see Eq. ( IA.33b . 

SCP-I 1 1 The values for g m ; n and g max for SCP-III are identical 
to those for SCP-II, and are given by Eqs. (IA.74b and ( IA.75t . 
The expression for the specific binary fraction for SCP-III is 
the same as in (the one for SCP-I). Note, however, that the pri- 
mary mass distribution /m, (Mi) is different due to the different 
limits g m in and g ma x- The right-hand panel in Fig. [7] shows an 
example of Sm, (Mi ) resulting from SCP-III. The term fu c (M\ ) 
vanishes for c < M\ < 2c, so that 2?m,(Mi) = 1 in this 
mass range. A discontinuity appears at Mi = 2c, beyond which 
Sm,(Mi) reaches its lowest point, after which it steadily rises 
until Mi = d. Beyond Mi = d, the specific binary fraction 
decreases again to Sm 1 (Mj ) = at M\ - 2d, for the same rea- 
sons as for SCP-I and SCP-II. The overall binary fraction S a n 
is always equal to S for SCP-III. 

A. 4. 3. Uniform mass ratio distribution and single 
power-law mass distribution 

In the expressions given in this section the term 1 - (c/d) y has 
often been absorbed in the normalisation constants (which are 
therefore not strictly consistent with the expressions above). 
The following is assumed for h q and f c (M c ) = f(x) (using x = 
M c for ease of notation): 

h q (q) =1 < q < 1 
f(x) = ax~ a 2c < x < 2d , 

where a = y(2c) 7 /(l - (c/d) 7 ) and y — a - 1, and a + 1. 

The SCP-I case Without restricting the primary mass range 
the expression for the mass ratio distribution is simply f q (q) - 
1. When the primary mass range is restricted the cases listed 
above have to be worked out. The corresponding expressions 
forfq(q) are listed in Table I A. fl Fig. lA.7l shows an example of 
what fq{q) looks like for all the SCP-I cases discussed above. 

The SCP-II case Without restrictions on the primary mass 
range the expression for the mass ratio distribution follows 



from Eq. (IA.6OI 1 and is: 

k 



/«(«) = 



1 - (c/d) 



2q_ 
1 + q 



(if 



(A.80) 



where the normalisation constant follows from: 
-1 

1 = I f q (q)dq 



f 



1 - (c/d) 



r (Znj i+ <? 



dq ■ 



■?mto) . (A.81) 



where g m ; n = (c/d)/(2 - c/d). The integral in this expression 
will be dealt with below. With a restricted primary mass range 
the expressions for f q are obtained for the three cases discussed 
above. These expressions are listed in Table IA.2I Again the 
term 1 - (c/d) y is absorbed in the normalisation constants k. 
All the normalisation constants for the SCP-II case contain the 
following integral: 

2q 



dq 



(A.82) 



For y + 0, -1, -2, -3, ■ ■ ■ the integral evaluates to an expres- 
sion involving a hypergeometric function: 



2y + VM(y+ l,y;y + 2;-g) 

y + 1 



(A.83) 



where the result was obtained by using the website 
integrals.wolfram.com. The value y = is not allowed 
as a — 1 was excluded. For y = -1 (a — 0) the integral is: 



Uq + \nq). 



(A.84) 



From the properties of the hypergeometric function (see PCP- 
III case) it follows that for y+2 = -n the expression above does 
not converge. This means that for y = -p, where p = 2, 3, 4, ■ ■ ■ 
the integral has to be evaluated separately. In this case we can 
write: 



= 2-'J(l + i) dq 

= 2-*jXo(?)<r^> 



which evaluates to: 

2-" 



h^Z(?)TVi 



(A.85) 



The SCP-III case Without restrictions on the primary mass 
range the expression for the mass ratio distribution is: 



f q (q)= [ (^—p\ax a dx, (A. 

Jc(M/q) \ X— LCI 



86) 



where the term for h' q follows from h' q (q) = 1/(1 - c/(x - c)). 
Note that the distribution is normalised. To bring out better the 
dependence on c/d one can also write (substituting z = x/d): 



y(2c/dy r 

m) ~ l-(c/d)r I 



c/d)(M/q) 



Z - (c/d) 
Z - (2c/d) 



z' a dz (A.87) 
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The integral over z: 



fl z-(c/d) \ 
J \z-(2c/d)j 



z- a dz 



can again be evaluated using the hypergeometric function for 
a + 0, 1, 2, 3, • ■ ■, and the result is: 

Ti 77 — L • (A.88) 

2(a- 1) 

The excluded values of a again follow from the convergence 
properties of the hypergeometric function. For a = the inte- 
gral evaluates to: 

z-(c/d)ln(z-2c/d). (A.89) 

The case a = 1 is excluded and for a = p = 2, 3,4, • • • the 
integral can be written as: 



1 c/d 



dz (A.90) 



z>'- 1 (z-2c/d) zP(z-2c/d) 
The solution can be found again from formula 2.117(4) in 



Gradshteyn & Ryzhik (2007): 



V 1 . (2cV~ P ]n ( z-2c/d \ _ 

£>r=l (j,-l-r)(2c/d) r z"- , - r T W ) \ I ) , A n-^ 

/ n \ (A.91) 

With restrictions on the primary mass range the integral dis- 
cussed above has to be evaluated for the various integration in- 
tervals corresponding to the cases discussed for SCP-II. The 
normalisation constants k can be obtained from a numerical in- 
tegration of the functions f g (q). For x\ = X2 the mass ratio dis- 
tribution becomes f q (q) = 1/(1— c/(xi— c)) = {x\—c)l(x\ — 2c). 
Figs. |A.9l and lA.10l show examples of what f q {q) looks like for 
all the SCP-II and SCP-III cases discussed above. Finally, we 
have listed the expressions for the binary fraction as a function 
of primary mass in Table IA. 31 
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SCP-I 



Mass ratio distribution 



Normalization constant 



Case 1 



Mi) 



< q < qo 

k[l-(2c/x 2 y(l + q)-v] q <q<l 



1-7 



(1 - q )(l - y) - V-y(2c/x 2 y + 2c/x 2 



Case 2 



f q (q) =k[l- (2c/x 2 y(l + q)-?] < q < 1 



1-7 



(\-y)-{2clx 2 y(V-y - 1) 



Case 3 



f q (q) 



k ! f ', , w " < q < q 2 



qi < q < 1 



l 



1 - (c/dy 



92 - 



(2c/x 2 y((i+q 2 y-y-i) 



1-7 



+ l-q 2 



Case 4 



< q < qo 

£[l -(|) r (l +<?)-''] q <q< qi 



k = 



1-7 



2'-y\{2clx x y - (2c/x 2 y] - y(2clx i - 2cjx 2 ) 



Case 5 



.4(!) 7 -(!) r ]( 1 + ^ 



Jfc = 



i-y 



9l(l-y)+[(2c-/ii))'-(2c/.V2F]2 1 -r+(2c/A-2)r-2c-/i 1 



Case 6a 



/,(<?) 



fe [ 1 -(|) r ( 1 + «)" r ] o< ?<?1 

fe [(f) 7 -(l) r ]( 1 + ^ 

*[(|) r (l+?) W -(§)1 ?2<?<1 



i-y 



q,(l-y)-(2c/x 2 )r[(l+q0 1 - y -l] 

tl + 

[(2 c /*,)r-(2 c /* 2 F][(l +<?2 )i-r-(i +1?1 )]-y] 

l_y 

(2 r /.v,)r(2 1 -Ml+< 72 ) 1 -'')-(c/<iF(l-< ?2 )(l- r ) 



Case 6b 



/ ? (<?) 



£[1 - (3 ) 7 ] qi<q<qi 
fe[(S) 7 ( 1+ ?) 1 "" y -(3) 7 ] 





<7 2 (I-y)-(2c-/A- 2 F[(l +<?2 )l-r-l 

i-y 



(l-(c/i)?X«l-m)(l-y) 

i-y 



+ 



(2cAv 1 )r(2i-r-(l+ 9 ,) 1 -r)-( c /rf)!'(l-, 1 )(l- r ) 



Case 7 



/,(?) = *(! + $)" 



1-7 
2 1 "? - 1 



Case 8 



f q (q) : 



4(l) r -(l) r ]( 1+ ^ y 0<?<?2 

4(l) y ( 1+ ^-(^)1 ^-<q<l 





[C2c/*i)r-<3<*a)i'] [ ( i +82 )i-r-i] 

hz 

(2c/.v 1 )y(2i->'-(i+ g2 ) 1 -y) 

i-y 

(r/rfFU-^Xl-y) 



Case 9 



/,(<?) 



fe [(l) r -(l) r ]( 1+ ^ r 0<q<q 2 
*[(|) r d+?) 1 - r -(§)1 <Z2<<7<<73 

93 < 9 < 1 



W 

[(2c-/A 1 )T-(2c-/A 2 ))'J[(l+ g2 )l-y-lJ 



(2c/.v,))'((l+ 93 )i-r-(l +92 )i-7) 

i-y 

(c/*y(g 3 -m)(i-r) 



Table A.l. The specific mass ratio distribution f q {q) for pairing function SCPT, resulting from a power-law generating mass 
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SCP-II 



Mass ratio distribution 



Normalization constant 



Case 1 



/,(«> 



o <q< J- 



Case 2a 



/«(<?) 



< q < ± 

i*ra r -(sn ***** 



< — 



(i + qi ) l -y - (i + ^)'" r ] + 
2^-(i + rf1-(5) r Ci-*) 



'c/x 2 

i-y 

(2f/.vi) r 

i-r 



Case 2b 



fM 



< o < 

AT 

4(ir-(i) r ]( i+ ^ ^«<» 
«im-m ***** 



( ^,)Mw j (1 + g2)1 - y _( 1 + ^^ 

£Z^[2^-(l + ^]-(§) r (l- g2 ) 



Case 3a 



0<^< 



4(^) r -(f) r ] *>*«<£ 
t**<* 3 

q 3 < q < 1 



i-r 

(2c/.vi)'' 

i-y 



(l+ ?2 )'-''-(l + ^) W 
(l + 93 ) 1 -''-(l + ^) 1 " r 



-(§)"(*-«.) 



Case 3b 



/,(«) 



0<q<-^ 

? 3 < 3 < 1 



= X/,r (^) <**- 

^[(^r-(^n+ 

SMrM [ ( i + ?2 )i-r _ (i + s.^ + 



Table A.2. The specific mass ratio distribution for pairing function SCP-II, resulting from a power-law generating mass 
distribution f m (m) and a uniform generating mass ratio distribution h q (q). See Appendix IA. 4. 3l for details. 
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Mass ratios for SCP-I a= 2.35, c=0.200, d= 2.000 



Mass ratios for SCP-I a= 2.35, c=0.200, d= 2.000 





0.7 



0.4 0.6 

q 

Mass ratios for SCP-I a= 0.30, e=0.200, d= 2.000 



1.0 



0.0 



0.2 



0.4 0.6 

q 



o.s 



1.0 




Mass ratios for SCP-I a= 0.30, c=0.200, d= 2.000 
-i ■ ■ i i ■ -yr ■ ! r- 



0.4 0.6 

q 

Mass ratios for SCP-I a=-2.50, c=0.200, d= 2.000 



1.0 




0.4 0.6 

q 

Mass ratios for SCP-I ce=-2.50, c=0.200, d= 2.000 



1.0 




Case l:x,= 0.10, x 2 = 0.35 

q„= 0.14, q,= 3.00, q,=10.43, q 3 =39.00 

Case 2: x,= 0.10, x 2 = 1.00 

q„=-0.60, q,= 3.00, qj= 3.00, q 3 =39.00 

Case 3: x,= 0.10, x,= 3.00 

q =-0.87, q,= 3.00, q 2 = 0.33, q 3 =39.00 

Case 4: x,= 0.25, x,= 0.35 

q = 0.14, q,= 0.60, q 2 =10.43, q 3 =15.00 

Case 5: x,= 0.25, x 2 = 1.00 

q„=-0.60, q,= 0.60, q 2 = 3.00, q 3 =15.00 



Case 6a: x,= 0.25, x 2 = 2.20 
q =-0.82, q,= 0.60, q 2 = 0.82, q,=15.00 
Case 6b: x,= 0.25, x,= 3.50 
q =-0.89, q,= 0.60, q 2 = 0.14, q 3 =15.00 
Case 7: x,= 0.50, x 2 = 1.80 
q„=-0.78, q,=-0.20, q 2 = 1.22, q 3 = 7.00 
Case 8: x,= 0.50, x,= 3.80 
q =-0.89, q,=-0.20, q 2 = 0.05, q^ 7.00 
Case 9: x,= 2.10, x 2 = 2.30 
q„=-0.83, q,=-0.81, q 2 = 0.74, q 3 = 0.90 



Fig.A.7. Left: Mass ratio distributions f q {q) for split-core pairing I for a uniform generating mass ratio distribution and core 
masses from single power-law mass distribution, where 2c < M c < 2d. The values of a are 2.35, 0.30, and -2.50. The black 
solid curve shows the complete mass ratio distribution (for all binaries in the population) for c = 0.2 and d = 2.0. The other 
curves show what happens to the observed f q (q) if the primary mass M\ is restricted to x\ < M\ < xi- The ten cases for SCP-I 
from section lAAll are listed in the legend. Right: Same, but with the vertical scale changed to bring out some of the details in 
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Mass ratios for SCP-II, single power-law IMF: ct= 2.35 




Mass ratios for SCP-II, single power-law IMF: cc= 0.30 
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Mass ratios for SCP-II, single power-law IMF: ot=-2.50 
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Mass ratios for SCP-III, single power-law IMF: a= 2.35 
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Fig. A.9. Mass ratio distributions f q (q) for split-core pairing II and III for a uniform generating mass ratio distribution and core 
masses from single power-law mass distribution for a = 2.35, 0.30, -2.50, where 2c < M c < 2d. The curves are shown for 
c/d = 0.01, 0.1, 0.3, 0.5. 
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Mass ratios for SCP-II a= 2.35, c=0.200, d= 2.000 



Mass ratios for SCP-III a= 2.35, c=0.200, d= 2.000 




Case 1: x 1= 0.25, x 2 = 1.50, c/x,= 0.80, c/x 2 = 0.13 

q =-0.73, q,= 0.60, q 2 = 1.67, q 3 =15.0O 

Case 2a: x,= 0.30, x 2 = 3.80, c/x,= 0.67, c/x 2 = 0.05 

q =-0.89, q,= 0.33, q 2 = 0.05, q 3 =12.33 

Case 2b: x,= 0.30, x 2 = 2.10, c/x,= 0.67, c/x 2 = 0.10 

q„=-0.81, q 1= 0.33, q 2 = 0.90, q 3 =12.33 

Case 3a: x,= 2.10, x 2 = 3.80, c/x,= 0.10, c/x 2 = 0.05 

q„=-0.89, q,=-0.81, q 2 = 0.05, q,= 0.90 

Case 3b: x,= 2.20, x,= 3.00, c/x,= 0.09, c/x,= 0.07 

q„=-0.87, q,=-0.82, q 2 = 0.33, q,= 0.82 



Fig. A.10. Mass ratio distributions f q (q) for split-core pairing II and III for a uniform generating mass ratio distribution and core 
masses from a single power-law mass distribution with a = 2.35, where 2c < M c < 2d. The black solid curve shows the complete 
mass ratio distribution (for all binaries in the population) for c = 0.2 and d = 2.0. The other curves show what happens to the 
observed f q (q) if the primary mass M\ is restricted to x\ < M\ < x-i. Results are shown for a = 2.35 (top), a = 0.3 (middle) and 
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Table A.3. The specific binary fraction Sm, (Mi) for the nine pairing functions described in this paper, resulting from a power- 
law mass distribution /m(M) and a uniform mass ratio distribution h q {q). See Appendices lA. U - lA~4l for a detailed description, and 
Fig.|7]for a visualisation. 



